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SECTION  1 


i 

INTSOOUCTION 

t 

I 

* 

Luir  cooling  of  neutral  particle  betas  provides  a  promising  j 

sechsmsa  for  reducing  the  divergence  of  the  relativistic  atoeic  J 

betas  presently  being  considered  for  inclusion  in  space-baaed  j 

weapons  systsas  The  nature  of  the  cooling  process  ie  such  that 
the  eaittance  of  the  resulting  beu  is  not  subject  to  the 
constraints  of  Liouvi lie’s  theorea,  as  are  magnetostatic  or 
electrostatic  focusing  systeas  In  static  focusing  systsas,  the 
eaittance  of  the  beaa,  which  is  proportional  to  the  product  of 
the  beaa  sue  and  the  beaa  divergence,  is  a  conserved  quantity 
Thus,  the  only  way  to  reduce  beaa  divergence  is  at  the  expense  of 
beaa  sise  Laser  cooling  permits  beaa  divergence  reduction 
without  increasing  beaa  size  Smaller  beaas  permit  smaller 
optics  in  the  accelerator,  thus  reducing  the  size,  weight,  and 
cost  of  eventual  weapons  systeas 

Interest  in  laser  cooling  technologies  has  also  recently  been 
expressed  in  the  related  area  of  atomic  trapping  Neutral 
particles  can  be  trapped  in  an  electromagnetic  field  and  confined 
there  for  extended  periods  of  time  Applications  of  such  trapped 
particles  include  precision  atomic  spectroscopy  and  containment 
and  study  of  anti -matter  Atomic  spectroscopy  has  DoD 
applications  in  precise  atomic  clocks  and  frequency  standards. 

Anti-matter  containment  is  important  in  concepts  currently  being 
considered  by  the  Air  Force  for  use  in  exotic  fuels  and  weapons.1 

Under  this  contract,  we  studied  three  areas  relevant  to  the 
cooling  and  trapping  of  neutral  atoms  in  electromagnetic  fields. 

The  first  of  these  areas  was  the  Optical  Kapi tza-Di rac  Effect 
(often  referred  to  as  the  Optical  Stern-Gerlach  Experiment)  which 
concerns  the  result  of  the  interaction  of  an  atomic  beam  with  a 
resonant  or  nearly  resonant  standing  wave.  The  onset  of  this 
effect  provides  one  of  the  fundamental  limits  in  one  of  our 


proposed  cooling  schemes  In  addition,  the  effect  provides  an 
interesting  aani f sstation  of  the  fundamental  interaction  between 
atoms  and  light. 

In  the  second  area  of  study,  we  measured  the  two- photon 
ionisation  cross  section  of  atomic  hydrogen  at  1216  A,  the  Lyman- 
alpha  (L^)  resonant  transition  from  the  ground  state  (Is)  to  the 
first  excited  state  (2p)  This  two  photon  cross  section  will 
limit  the  maximum  cooling  rates  possible  for  a  beam  of  hydrogen 
atoms  If  the  laser  beam  is  too  intense,  a  significant 
proportion  of  the  neutral  atoms  will  be  ionised  during  the 
cooling  process 

Our  third  area  of  study  under  this  contract  the  development 
of  the  conce  pt  of  »  Pulse  Cooling  Under  appropriate 
circumstances,  w-pulse  provides  a  means  for  faster  cooling  than 
is  possible  with  the  conventional  steady-state  radiation  method 

Section  2  of  the  report  is  a  review  of  the  physics  of  the 
interaction  of  laser  radiation  with  neutral  atoms,  as  applied  to 
the  studies  done  in  this  program  The  section  emphasis  is  on  how 
radiation  forces  can  be  used  to  cool  and  trap  neutral  atoms.  The 
following  three  sections  discuss  our  research  in  the  areas  of  the 
Optical  Kapi tza-Di rac  Effect  (OKDE) ,  hydrogen  resonant  two-photon 
ionization,  and  w-Pulse  Cooling.  Section  6  summarizes  our 
research  results  and  puts  them  in  the  perspective  of  current 
research  efforts  in  laser  cooling  and  trapping.  Implications  of 
laser  cooling  and  trapping  for  current  and  future  DoD  projects  is 
mentioned  and  we  also  make  some  suggestions  for  further  research 
in  this  area.  The  final  section  is  a  list  of  publications 
stemming  in  whole  or  in  part  from  this  contract.  For  the 
reader’s  convenience,  we  have  included  the  full  text  of  each 
publication  - 


SECTION  2 


REVIEW  or  TO  P1YSICS  OF  ATOMIC  BEAM-LASER  INTERACTION 

When  an  atoa  interacts  with  an  external  electromagnetic 
field,  the  center  of  mass  motion  of  the  atom  is  modified  by  the 
recoil  of  the  atom  when  absorbing  or  emitting  photons  Under 
appropriate  conditions,  the  collective  atomic  recoil  serves  to 
reduce  the  momentum  spread  of  an  atomic  sample  (laser  cooling)  or 
to  produce  a  diffraction  pattern  from  a  beam  of  atoms  (the 
Optical  Kapi txa-Di rac  Effect) 

In  this  section,  we  will  review  some  physics  necessary  for 
understanding  the  most  prominent  features  of  these  two  effects 
In  addi  t  ion,  we  will  present  several  techniques  that  could  be 
used  to  cool  an  atomic  sample  The  theory  of  the  Optical 
Kapi tza-Di rac  Effect  will  be  developed  in  Section  3  Details  of 
the  specific  theoretical  and  experimental  developments  achieved 
with  the  support  of  this  contract  will  be  discussed  in 
Sections  3,  4,  and  5 

The  Hamiltonian  for  the  interaction  of  an  atom  with  a  laser 
field  can  be  written  in  the  form 

H  --  H0*  Hp*  Hcm  (2  1) 

H0  is  the  Hamiltonian  for  the  internal  atomic  degrees  of  freedom 
For  the  purposes  of  the  following  discussion,  we  will  assume  that 
the  unperturbed  atom  can  be  accurately  modeled  as  consisting  of 
two  levels,  a  ground  state,  | a> ,  and  an  excited  state  |b>  (see 
Figure  1).  Hj  represents  the  coupling  between  the 
electromagnetic  fields  and  the  atom.  In  the  cases  we  will 
consider  here,  the  dominant  contribution  to  Hj  is  given  by 
Hi  =  ~n* E,  where  p  is  the  atomic  dipole  moment  and  B  is  the 
electric  field.  Hp  is  the  free-field  Hamiltonian  for  the 
electromagnetic  field.  Hc H  represents  the  center  of  mass  motion 
of  the  atom,  where  HCM  =  p2/2m. 


In  the  low  intensity  limit,  the  eigenstates  of  the  atomic  and 
electromagnetic  fields  are  decoupled,  and  the  internal  and 

external  atomic  states  are  independent  of  the  intensity  of  the 

electromagnetic  field.  In  the  high  intensity  limit,  the 
electromagnetic  field  serves  to  mix  the  unperturbed  atomic 
eigenstates.  This  mixing  can  be  most  easily  handled  by  replacing 
the  "pure"  atomic  states  with  "dressed  states"  (eigenstates  of 
the  atom  plus  laser  field  system),  which  will  be  denoted  as  |+,n> 
or  |-,n>,  where  n  is  the  photon  number  of  the  state. 

The  dressed  states  occur  in  doublets,  each  state  consisting 

of  a  mixture  of  the  two  states  |a,n+l>  and  |b,n>.  Let  6  be  the 
detuning — the  difference  between  the  laser  frequency  wL  and  the 
transition  energy  Eb-E, .  We  will  denote  the  Rabi  flopping 
frequency  as  u/R  ,  where 


fa/p  =  n»E  /h 


(2.2) 


The  quantity  fi  is  the  dipole  interaction  matrix  element  e<b|r|a> 
and  E  is  the  electric  field  vector.  If  the  phase  of  the  dipole 
interaction  matrix  element  is  chosen  to  be  real,  the  amount  of 
mixing  of  the  two  states  is  governed  by  the  intensity  dependant 
expressions  4 


+  ,n>  =  cos  (6)  l  a  ,  n  ■»  1  >  -*  sin(d)  b,n> 
,n>  -  sin(0)!a,n-*l>  ♦  cos  (fi)  b  ,  n  >  , 


(2  3) 


where  the  mixing  angle  6  is  defined  by  tan(0)  (fi  6 ) ■ w*  and  fl  is 
defined  to  be  *  The  energies  of  the  two  eigenstates 


(  n  *  1  i  u/, 


■  n  ♦  1  J  w, 


A  number  of  schemes  have  been  used  to  cool  a  gas  by  taking 
advantage  of  the  exchange  of  momenta  between  photons  and  the 
constituent  atoms  of  the  gas.  Four  general  schemes  have  been 
proposed,  which  we  will  designate  as  collisional,  spontaneous, 
stimulated,  and  T-pulse.  Realizations  of  each  of  these  schemes 
will  be  discussed  below. 

1.  Collisional  Cooling 

An  example  of  a  collisional  cooling  scheme  is  shown  in 
Figure  2.  A  laser  is  tuned  to  the  red  side  of  an  atomic 
resonance,  exciting  the  atom  from  | a>  to  | b> .  The  atom 
spontaneously  decays  from  | b>  to  | a> ,  losing  an  energy  Eb-E, , 
slightly  more  than  the  energy  of  the  exciting  laser.  The  extra 
energy  is  lost  in  collisions  with  other  atoms-  in  the  gas.  This 
technique  is  the  least  efficient  of  the  laser  cooling  prospects 
we  will  discuss.  The  change  of  momentum  per  absorbed  photon  is 
typically  only  a  fraction  of  the  photon  momentum. 

2.  Spontaneous  Cooling 

A  much  more  efficient  technique  for  cooling  atoms  is  what  we 
call  spontaneous  cooling  A  group  at  the  National  Bureau  of 
Standards  was  able  to  use  tnis  technique  to  cool  a  thermal  beam 
of  sodium  atoms  down  to  a  temperature  of  100  mK  ‘  Figure  3  :s  an 
illustration  of  the  technique  we  originally  envisioned  foi 
transverse  coo  ling  of  a  relativist  it  a t  om i '  t  ■  e  am  t  >  spontaneous 
coo  1  l  ng  In  this  technique,  two  »  as  e  r  he  ams  are  .  r.  i  it-ti  t  or, 
either  side  of  the  a  torn  l  <  team  to  re  .  -><  ,  erj  T  r.  *  frequency  i  ♦ 
the  iight  seen  in  the  res’  fr  aae  of  e  a  r  a  t  oil  *  .  .  .  t  e  [>opp  .  e  r 
shifted  to  a  new  frequen  >  e  wh  i  «  h  'legends  •  he  »  .  ...  i  ♦  v 

the  atom  If  0  is  the  s  [  A  ’he  a’  >o.  -  '  i  r.  .  *  s  '  'he  a  pee  'l  f 

light,  and  7  1  0 ‘  *  1  ’he  frequen' v  .  is  g  .  v  e:  r \ 


71  1  0  os  t> 


.  V.  h  ,  '  *  *h  i  *  . 
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Figure  3.  Spontaneous  cooling  scheme.  Atoms  moving  towards 

laser  beam  1  are  resonant  with  it,  but  non-resonant 
with  laser  beam  2.  They  are  "pushed”  away  from 
laser  1  by  the  recoil  of  the  absorbed  photon  from 
laser  1 .  Since  the  atom  scatters  the  absorbed  photon 
into  4?r,  the  net  effect  of  the  fluorescence  is  0. 
Atoms  moving  toward  laser  beam  2  are  likewise  moved 
toward  a  central  trajectory. 


As  indicated  in  Figure  3,  9  is  the  angle  between  the  laser  beam 
and  the  atomic  velocity  vector.  In  the  non-relativistic  limit  of 
velocities  much  slower  than  the  speed  of  light,  this  expression 
reduces  to 

» 

v  -v  -  -k*v/2ir  ,  (2.6) 

where  k  is  the  wave  vector  of  the  laser  beam. 

Spontaneous  cooling  techniques  take  advantage  of  this 
inhomogeneous  Doppler  shift  to  cool  the  atoms  in  the  beam.  For 
simplicity,  imagine  a  beam  in  which  all  of  the  atoms  have  a 
uniform  speed,  but  slightly  different  trajectories.  The  cos(0) 
term  of  Eq.  (2.5)  will  cause  the  resonant  frequency  of  each  atom 
to  have  an  angle-dependent  frequency  shift.  Thus  the  laser  will 
only  be  resonant  with  atoms  within  a  narrow  angular  range.  These 
atoms  will  pick  up  one  unit  of  photon  momentum,  hk,  as  they 
absorb  a  photon  from  the  laser  beam.  An  equal  amount  of  momentum 
is  lost  when  the  photon  is  reradiated  as  the  atom  relaxes  to  the 
ground  state.  However,  the  photons  absorbed  from  the  laser  are 
all  in  the  same  direction  and  their  effect  is  cumulative.  The 
photons  radiated  by  the  excited  atoms  are  radiated  into  a  4* 
solid  angle.  The  effects  of  the  radiated  photons  average  out  to 
zero . 

In  the  setup  illustrated  in  Figure  3,  atoms  moving  away  from 
the  central  trajectories  of  the  atomic  beam  will  be  "pushed"  away 
from  the  lasers  they  are  moving  toward.  Atoms  moving  on  the 
central  trajectory  or  away  from  a  particular  cooling  laser  will 
not  be  non-resonant  with  the  laser  light  and  therefore 
unaffected.  As  the  cooling  process  continues,  the  distribution 
of  atomic  resonance  frequencies  will  shift  toward  the  frequency 
of  the  atoms  with  the  central  trajectory.  To  further  cool  the 
distribution,  the  laser  frequency  needs  to  be  shifted  to  remain 
resonant  with  atoms  on  the  periphery  of  the  velocity 
distributions.  Alternatively,  an  external  electric  or  magnetic 


field  or  both,  could  be  applied  to  slightly  shift  the  atomic 
resonance  via  the  Stark  or  Zeeman  effects  as  the  distribution  is 
cooled . 

The  fact  that  all  real  beams  have  a  non-zero  spread  in  atomic 
velocities  slightly  complicates  the  process  of  spontaneous 
cooling,  because  the  frequency  shift  of  the  atoms  depends  on 
their  speeds  as  well  as  on  their  trajectories.  For  example, 
consider  a  beam  of  250-MeV  hydrogen  atoms  with  a  gaussian 
momentum  spread  ap=pxlO*4,  which  is  being  cooled  by  means  of  the 
La  transition  at  1216  A.  Assume  that  the  angular  divergence  of 
the  beam  is  one  mrad .  If  the  cooling  laser  is  oriented  75*  with 
respect  to  the  atomic  beam,  the  velocity  variations  in  the  beam 
result  in  a  variation  of  60  GHz  in  the  required  resonant  laser 
frequency.  The  trajectory  variations  in  the  beam  will  require 
only  a  1.6-GHz  frequency  spread.  In  short,  there  is  a  problem  in 
that  the  frequency  shifts  caused  by  the  divergence  could  be 
masked  by  those  caused  by  longitudinal  momentum  spread. 

Fortunately,  this  problem  can  be  solved  by  an  appropriate 
orientation  of  the  cooling  lasers.  If  the  angle  between  the 
beams  is  chosen  to  be  the  so  called  "magic  angle"  Qm,  given  by 

6n  =  arcos  (/?)  ,  (2.7) 

the  effects  of  speed  variations  on  the  frequency  spread  in  the 
atoms  will  be  minimized.  In  the  example  above,  the  velocity 
variations  in  the  beam  will  then  result  in  a  variation  of  only 
8.6  MHz  in  the  laser  frequency.  The  trajectory  variations  will 
demand  a  variation  of  2.5  GHz.  With  a  non-relativistic  beam,  the 
magic  angle  reduces  to  00*,  as  would  be  expected  from  Eq.  (2.7). 

2.  Stimulated  Cooling 

A  disadvantage  of  the  spontaneous  cooling  mechanism  outlined 
above  is  that  the  minimum  cooling  time  is  limited  to 


'l 


r  =  27fik/v 


perp 


(2.8) 


where  7  is  the  spontaneous  decay  lifetime,  k  is  the  photon  wave 
vector,  and  vp#rp  is  the  maximum  transverse  velocity  component  of 
the  atomic  beam.  With  an  intense  laser  beam,  it  is  possible  for 
the  stimulated  decay  rate  of  the  atoms  to  greatly  exceed  the 
spontaneous  decay  rate.  Aspect  et  al . 3  have  demonstrated  such  a 
cooling  scheme4,  utilizing  stimulated  emission  that  reduced  the 
cooling  time  by  an  order  of  magnitude  over  that  achieved  by  the 
spontaneous  process.  Their  non-relativistic  approach  consists  of 
passing  the  atoms  through  a  standing  wave  that  is  tuned  to  the 
high  energy  side  of  an  atomic  resonance  frequency.  The  standing 
wave  provides  a  periodic  spatial  intensity  distribution  that 
induces  a  spatial  periodicity  in  the  energies  of  the  atomic 
states,  as  illustrated  in  Figure  4.  The  admixture  of  unperturbed 
atomic  states  in  the  dressed  |+>  and  |->  states  also  shows  a 
position  dependence. 

Consider  an  atom  with  a  velocity  that  has  a  small  transverse 
component  along  the  periodicity  direction  of  a  standing  wave.  If 
the  atom  begins  in  the  |  +  ,n>  state  at  a  node  of  the  standing 
wave,  cos(0)=l  and  the  atom  is  in  a  pure  |a,n+l>  state,  and  at 
the  node  of  a  transverse  potential  well.  As  the  atom 
adiabatically  moves  to  an  anti-node  in  the  distribution,  some  of 
the  |b,n>  state  is  mixed  into  the  configuration.  Spontaneous 
decay  to  a  |+,n-l>  or  a  |-,n-l>  becomes  possible.  Upon  decay  to 
a  ]-,n-l>  state,  the  atom  will  again  find  itself  at  the  bottom  of 
a  transverse  potential  well  that  it  will  climb  again.  The  energy 
required  to  climb  the  potential  well  creates  a  drag  force,  which 
serves  to  damp  out  the  transverse  velocity  component  of  the  beam. 

This  technique  shows  promise  for  greatly  enhanced  cooling 
rates  without  saturation.  Further  research  should  be  pursued  how 
this  method  could  be  applied  to  cooling  a  relativistic  beam. 
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Figure  4.  Position-dependence  of  "dressed"  atom  energies  and 
states  in  a  standing  wave. 
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8.  w-Pulse  Cooling 


Another  way  to  improve  on  the  cooling  rate  limit  of 
spontaneous  cooling  is  to  stimulate  radiative  decay  rather  than 
wait  for  the  excited  state  to  decay  spontaneously.  By 
appropriate  selection  of  the  duration  and  intensity  of  the 
cooling  pulses,  it  is  possible  to  obtain  a  unit  probability  of  a 
transition  from  the  ground  state  to  the  excited  state,  or  vice 
versa.  Such  a  pulse  is  called  a  x-pulse .  Cooling  an  atomic  beam 
with  this  technique  is  illustrated  in  Figure  5.  The  atom  is 
originally  in  the  ground  state  when  it  is  hit  with  the  first 
K  pulse.  This  pulse  gives  the  atom  a  unit  probability  of  being 
in  the  excited  state,  and  transfers  one  unit  of  photon  momentum 
to  the  atom.  A  second  ir-pulse  is  then  applied  from  the  opposite 
side  of  the  beam  to  stimulate  the  atom  back  down  to  the  ground 
state.  The  atomic  recoil  from  the  second  pulse  is  in  the  same 
direction  as  the  momentum  picked  up  from  the  first  pulse.  If  the 
Doppler  effect  is  used  to  act  selectively  on  only  those  atoms 
moving  towards  the  source  of  the  first  pulse,  this  process  is 
very  similar  to  the  spontaneous  cooling  scheme.  The  w-nulse 
cooling  rate  will  be  given  by  the  stimulated  decay  rate,  which  is 
much  faster  than  the  limiting  spontaneous  decay  rate  of  the 
spontaneous  cooling  process. 

Further  details  of  this  process  can  be  found  in  Section  5  and 
in  the  Palmer  and  Lam  paper  included  in  the  Publications  section. 


13616-9 


-  lb  >  »  lb  >  -  lb  > 

TWO- LEVEL  PI-PULSE  PI-PULSE  _ 

ATOM  - ► 


|6> 


l»> 


P1  '  0 


Pi  *  hk 


PI -PULSE 

lPl 


2hk 


PI -PULSE 


Figure  5.  ir-pulse  cooling  scheme.  Atoms  are  hit  with  a  w-pulse 
while  in  the  ground  state,  giving  them  one  unit  of 
photon  momentum.  The  atoms  are  then  stimulated  back 
to  the  ground  state  from  the  opposite  side,  giving 
them  another  unit  of  photon  momentum  in  the  same 
di rection . 
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SECTION  8 


OPTICAL  KAPITZA-DIRAC  EFFECT 


A.  INTRODUCTION 

The  Optical  Kapitza-Dirac  Effect  occurs  when  a  beam  of  atoms 
passes  through  a  standing  wave  region  of  resonant  or  near 
resonant  radiation.  A  physical  picture  of  many  of  the  aspects  of 
this  effect  can  be  obtained  by  considering  the  energy  levels  of 
the  dressed  atomic  states  in  a  spatially  modulated  laser  field. 

If  the  atoms  enter  the  interaction  region  adiabatically  from  the 
ground  state,  they  will  stay  in  the  |+>  state  during  the 
interaction  (in  the  absence  of  spontaneous  decay) .  The  gradient 
of  the  field  intensity  results  in  a  position-dependent  force 
being  exerted  on  each  atom.  An  initially  well-collimated  beam 
will  acquire  a  transverse  velocity  distribution  as  illustrated  in 
the  smooth  curve  of  Figure  6.  A  full  quantum-mechanical 
treatment  of  the  effect  results  in  a  curve  consisting  of  discrete 
lines,  as  shown  in  the  second  curve  of  the  figure.  The  discrete 
spectrum  results  from  the  fact  that  momentum  transfer  to  the  atom 
is  always  in  the  direction  of  the  laser  beam,  and  occurs  in 
discrete  units  of  the  photon  momentum. 

A  semi-classical  calculation  of  this  effect  was  done  first  by 
Cook  and  Bernhardt.  Later  refinements  by  Bernhardt  and  Shore,® 
Arimondo  et  al . , 7  and  by  Tanguy  et  al . 8 ’ 9  involved  quantizing  the 
radiation  field,  detuning  from  resonance,  considering  the  effect 
of  spontaneous  emission,  developing  a  Gaussian  laser  beam 
profile,  and  misaligning  the  laser  and  atomic  beams. 

Concurrent  work  on  this  effect  has  been  done  by  Pritchard’s 
group  at  M.I.T.,  in  which  they  obtained  similar  results  to  those 
reported  here.1® 


B.  THEORETICAL  STUDIES 


The  current  understanding  of  spontaneous  radiative  cooling  is 
that  the  optimum  intensity  for  cooling  is  equal  to  the  saturation 
intensity  of  the  atomic  transition.  This  conclusion  raises  the 
important  question  of  how  is  the  atom’s  translational  motion 
perturbed  by  such  a  strong  laser  field.  Past  studies  of  atom- 
laser  interactions  point  to  the  fact  that  in  the  regime  where  the 
laser  intensity  is  equal  to  or  greater  than  the  saturation 
intensity  of  the  atom,  the  atomic  structure  is  significantly 
modified,  and  it  is  not  possible  to  treat  the  atom  and  the 
radiation  field  as  separate  entities.  Therefore,  we  must 
approach  the  problem  from  the  quantum  electrodynamics  (QED)  point 
of  view.  The  essence  of  this  approach  is  solving  the  entire 
problem  of  the  atom  plus  field  exactly  (i.e.,  finding  the 
eigenstates  of  the  total  Hamiltonian.  We  then  use  the  exact 
eigenstates  to  construct  experimental  observables.  For  atomic 
motion,  the  observable  is  the  momentum  distribution  of  the  atoms 
in  the  presence  of  the  radiation  field. 

This  technique,  called  the  dressed  atom  approach,  was 
pioneered  by  Cohen-Tannoud j i  in  his  studies  of  resonance 
fluorescence  in  intense  fields.11  The  term  "dressed”  atom  has  an 
intuitive  appeal  because  the  eigenstate  of  the  total  Hamiltonian 
contains  the  photon  number  states,  which  can  be  interpreted  as 
the  "cloth"  for  the  "bare"  atomic  states.  In  QED  terminology, 
the  physical  properties  of  the  atom  have  been  renormalized  by  the 
presence  of  the  photon  field. 

This  section  develops  the  dressed  atom  solutions  in  several 
stages.  First,  we  will  discuss  the  bare  atom  plus  photon  states. 
Then  we  will  introduce  the  concept  of  electric  dipole  coupling 
for  generating  the  dressed  states  for  a  unidirectional  photon 
field.  We  will  then  generalize  the  dressed  states  to 
bidirectional  photon  fields  (i.e.,  the  standing  wave  problem). 


Consider  the  two-level  system,  shown  in  Figure  1,  where  I a> 
and  I b>  are  the  ground  and  excited  states  respectively.  The 
Hamiltonian  of  the  bare  atom  plus  photon  field  is  given  by 

H0  =  p2/2M  +  E  hwa  ba*ba  +  E  fcu/k  (ak*ak  +  1/2)  ,  (3.1) 

where  p2/2M  is  the  kinetic  energy  of  the  center  of  mass  of  the 
atom,  h u)a  is  the  energy  of  the  atomic  state  la>,  and  ba  and  ba* 
are  the  respective  lowering  and  raising  operators  respectively. 

In  terms  of  the  Dirac  Ket  and  Bra  notation,  ba*ba  =  la><al .  The 
term  hwk  denotes  the  energy  of  the  photon  having  wavevector  k;  ak 
and  ak *  are  the  annihilation  and  creation  operators, 
respectively,  for  the  photon  field.  The  sum  over  k  is  over  all 
values  of  k  and  polarization  states  of  the  photon  field.  The 
eigenstates  of  H0  are  given  by 

lp,a,nk>  =  lp>  x  la>  x  lnk>  ,  (3.2) 

whe*-e  I  p>  is  the  eigenstate  of  p2/2M,  and  lnk>  is  the  eigenstate 
of  lnk>.  The  product  sign  indicates  all  possible  product 
combinations.  The  direct  products  of  "bare"  states  in  Eq.  (3.2) 
are  a  result  of  the  fact  that  the  Hamiltonian  does  not  contain 
the  coupling  between  the  atomic  and  photon  fields. 

To  generate  the  dressed  states,  we  start  with  the  Hamiltonian 
for  an  atom  interacting  with  a  single-mode  photon  field, 

Hi  =  p2/2M  +  E  hwa  ba*ba  +  V(t)  ,  (3.3) 

where  the  electric  dipole  coupling  to  the  radiation  field  is 
given  by 

V(t)  =  (1/2)  n»E0  bb’b,  a  exp(  ik»R  -  u/t) 

-(1/2)  M  Eo  b,*bb  a*  exp(-ik*R  4  wt)  ,  (3.4) 


:4 


and  where  E0  =  i (Rw/2e 0 v) 1 '  2  is  the  amplitude  of  the  photon 
field.  In  the  derivation  of  Eq.  (3.4)  we  have  invoked  the 
rotating  wave  approximation.  It  should  be  noted  that  the 
Hamiltonian  Hx  describes  the  atomic  system  only.  The  total 
Hamiltonian  is  generated  by  means  of  a  unitary  transformation, 
which  does  not  change  the  physical  description  of  the 
observables . 

The  Hamiltonian  Hj  satisfies  the  Schroedinger  equation 

ihd  l*>/8t  =  Ht  l*>  ,  (3.5) 

which  is  time  dependent  since  V(t)  depends  on  time  through  the 
phase  wt.  In  the  Schroedinger  picture,  the  eigenstates  of  an 
operator  must  be  time  independent.  Hence  it  is  desirable  to  work 
in  a  reference  frame  where  V  is  time  independent.  Let  us 
introduce  a  unitary  transformation,  U  =exp(iHft),  where 
Hf  =  fiu/(a*a  +  1/2)  ,  so  that 

l+>  =  U  l$>  ,  (3.6) 

where  I4>>  is  the  state  vector  in  the  new  frame  of  reference. 
Carrying  out  the  transformation  for  the  Schroedinger  equation 
[Eq.  (3.5)],  we  find  that  l$>  satisfies 

ih  ai*>/3t  =  H2  l*>  ,  (3.6) 

wi  th 

H2  =  H*e  -  (l/2)^iE0  b  b  *  b  ,  a  exp(ik«R) 

-(1/2)mE0  b,*bb  a’  exp ( -ik»R)  ,  (3.7) 

where  the  superscript  s  denotes  a  single  mode  radiation  field  in 
Eq.  (3.1)  . 

The  eigenstates  of  H2  can  be  found  by  solving  the  eigenvalue 
problem , 


where 


I  ♦>  =  exp(-iAt)  !±> 


(3.9) 


They  are  given  by 

!+,n>  =  sin# I p ’ ,b,n-l>  +  cos#lp,a,n>  (3.10) 

l-,n>  =-cos# I p ’ , b , n-l>  +  sin#lp,a,n>  ,  (3.11) 

with 

sin#  =  <b!H2la>/N  (3.12) 

cos#  =  (E  -  A,)/N  (3.13) 

N  =  {  (E  -  A.)2  +  I <b I H2 I a> I  2  }1/2  (3.14) 

<blH2la>  =  -(l/2)/iEc  (n)1/2  <p  ’  i  exp  ( i  k«  R)  I  p>  (3.15) 

E  -  (p ’ ) 2 /2M  ♦  hwb  *  hw(n  -  1/2)  (3.16) 

A.  =  (E-*G)  /2  *{  (hfi)2  *  4C  }  1  /  2  /  2  (3  17) 

G  =  p2/2M  +  hw,  +  hw ( n  ♦  1 / 2)  (3  18) 

C  -  (1/4)I/*E0I2  n  l<p  lexp(ik*R)!p>l2  (3  19) 

<5  =  {  (p  '  )  2  p2  >,  2MFi  •>  A  (3  20) 


with  p 


p  ♦  hk  by  conservation  of  1  inear  momentum 


Figure  7  shows  the  energy  levels  of  the  bare  atom  plus  field. 
Figure  8  shows  the  energy  levels  of  the  dressed  atom.  The  energy 
separation  of  the  dressed  states,  l+>  and  l->,  is  given  by 
Rfl  =  R (62+4C) 1 ^ 2 ,  the  generalized  Rabi  frequency.  It  should  be 
noted  too  that  the  atom  can  absorb  from,  and  emit  to,  the  photon 
field  many  times.  Hence  the  dressed  atom  consists  of  an  infinite 
ladder  of  pair  states  as  shown  in  Figure  9.  The  dressed  states 
[Eqs .  (3.10)  and  (3.11)]  are  the  main  results  of  this  section. 

We  can  write  the  eigenstate  equation  for  a  bidirectional 
photon  field  analogously  to  the  calculations  for  a  dressed  atom 
in  a  unidirectional  photon  field  as 

H3  I D ; >  =  A;  I D , >  ,  (3.21) 

where  the  Hamiltonian  H3  is  given  by 

H3  =  H0  *  ‘  -  -  (l/2)/xE0  bb'ba  { a  k e  1  k  *  R  +  a_ke-lk*R} 

“  (l/2)/iE0  b,-bb  {ak  *  e  '  k  *  R  +  a.k*e'ik*R}  ,  (3.22) 

and  where  the  subscripts  k  and  -k  denote  the  bidirectional 
wavevectors.  The  Hamiltonian  H0‘ ' * •  is  given  by 

H0*  *  =  p 2 /2M  +  T.  fiwQ  ba*bQ  +  Rw  (ak*ak  +  1/2) 

♦  Rw  (a_k’a_k  +  1/2)  (3.23) 

A  direct  computation  of  Eq .  (3.21)  gives  the  following  expression 

for  the  bidirectional  dressed  eigenstates: 

1 D ,  , n  k , n  _  k  >  =  sin^i  sin0;  lp',b,nk,n.k-l> 

+  cosf  sin0;  Ip  ,b,nk-l,n.k> 

+  COS e ■  I p , a , n  k ,n.k> 


(3.24) 


where 


cos 6 ; 

=  1/N-, 

(3.25) 

COS 6  ; 

sin#; 

=  <blH,  *la>/<  (F-A  •,  )  N ; 

> 

(3.26) 

Bind-, 

s in#  , 

=  <b 1 Hj " 1 a>/{  (E-A , )N; 

> 

(3.27) 

N 

=  {  1  +  1 <blHs ' 1 a> 1 2/(F-A ; ) 2 

+  1 <b 1 Hj  *  1 a> 1 2/ (E-A ; ) 2  > 

(3 . 28) 

E 

=  (p*)2/2M  +  hwb  +  hw 

(nk  +  n.k) 

(3.29) 

F 

=  (p’)2/2M  +  hwb  +  R w 

(“k  +  n.k) 

(3.30) 

G 

=  p2/2M  +  +  Rw  (nk 

♦  n.k) 

(3.31) 

<b  1  H 

l3’la> 

=  (1/2)ME0  bb*b#{(nk) 1 

1  2  <p ’ 1 e ' k  *  R 1 p> 

-*■  (n_  k )  1  /  2  <p ’ 1 e' ' k  * 

Rlp»  • 

(3.32) 

Here  the  expression  with  the  superscript  "  is  obtained  from  those 
with  superscript  ’  by  the  interchange  of  the  corresponding 
momentum.  The  dressed  state  energies  A  •,  are  roots  of  the  cubic 
equation , 

(E-A) (F-A) (G-A)  -  I <b I H j ’ I a> I  2  (E-A) 

-  I  < b I H 3 " I a> I  2  (F-A)  =  0  (3.33) 

Figure  10  shows  the  eigenstates  of  the  bare  atom  plus  the 
bidirectional  field.  The  frequencies  6  and  6*  have  expressions 
identical  to  those  defined  for  6.  The  superscripts  ’  and  "  stand 
for  the  momentum  p  =  p  +  hk  and  p*  =  p  -  hk ,  respectively. 

Figure  11  shows  a  set  of  triplet  dressed  states.  The  frequency 
separation  can  be  found  by  solving  the  cubic  Eq.  (3.33). 


Eqs .  (3.24)  constitute  the  main  result  of  this  section.  A  word 
of  caution  is  necessary  about  the  regime  of  validity  of  our 
solutions.  The  dressed  atom  solutions  are  valid  provided  that 
the  laser  field  is  tuned  off  resonance  by  more  than  a  natural 
linewidth.  Otherwise,  there  will  be  interference  effects  arising 
from  the  standing  wave,  leading  to  a  far  more  complicated 
solution  of  the  problem. 

In  summary,  we  have  presented  an  exact  description  of  the 
dynamics  of  two-level  atoms  in  the  presence  of  photon  fields. 

The  generalization  of  the  stationary  "dressed"  atom  approach  to 
moving  atoms  shows  that  the  dressed  energies  and  states  depend  on 
the  linear  momenta  of  the  atoms  and  photons.  For  a 
unidirectional  photon  field,  the  dressed  states  are  given  as  a 
linear  coherent  superposition  of  two  bare  atom  plus  photon  field 
states:  the  ground  state  of  the  atom  moving  with  momentum  p  and 
having  n  photons  present,  and  the  excited  state  of  the  atom 
moving  with  momentum  p’  =  p  +  Rk  and  having  n  -  1  photons 
present.  The  frequency  separation  between  the  two  dressed  states 
is  given  by  the  momentum-dependent  Rabi  frequency. 

In  bidirectional  photon  fields,  the  dressed  states  are  linear 
superpositions  of  the  three  bare  atom  plus  bidirectional  photon 
field  states:  the  ground  state  of  the  atom  moving  with  momentum 
p  and  having  nk  and  n_k  photons  present;  the  excited  state  of  the 
atom  moving  with  momentum  p’  =  p  +  hk  and  having  nk  -  1  and  n_k 
photons  present;  and  the  excited  state  of  the  atom  moving  with 
momentum  p  =  p  -  hk  and  having  nk  and 

n.k  -  1  photons  present.  The  frequency  separation  among  the 
three  dressed  states  is  determined  by  the  solution  of  Eq.  (3.33). 

These  solutions  constitute  the  starting  point  for  the 
calculation  of  the  momentum  distribution  of  the  atoms  in  the 
presence  of  a  resonant  photon  field.  An  inspection  of  the 
dressed  state  coefficients  shows  that  for  a  bidirectional  field, 


an  initially  symmetric  single  peak  momentum  distribution  function 
yields  a  symmetric  splitting  around  the  initial  maximum,  giving 
rise  to  the  Optical  Kapitza-Dirac  Effect.  Hence  the  dressed  atom 
picture  appears  to  provide  a  self-consistent  approach  toward  the 
understanding  of  the  effects  of  radiation  forces  on  atoms. 

C.  IXPBtHfBNTAL  STUDIES 

We  experimentally  tested  the  predictions  of  the  Optical 
Kapitza-Dirac  Effect  (OKDE)  theory  using  the  apparatus  whose 
block  diagram  is  shown  in  Figure  12.  Sodium  atoms  originate  in 
an  oven  and  are  collimated  by  the  snout  and  a  pair  of  slits.  The 
atoms  then  pass  through  a  fluorescence  monitor  that  is  used  to 
monitor  the  beam  density  and  the  laser  frequency.  In  the 
succeeding  interaction  region,  the  atoms  are  prepared  so  that  the 
only  atoms  present  in  the  F=2  ground  state  lie  within  a  very 
narrow  velocity  group  These  atoms  are  then  pumped  into  the  mf=2 
state  to  experimentally  prepare  a  "two  state"  system  The  atoms 
next  enter  the  standing  wave  region  where  the  OKDE  interaction 
takes  place.  The  spatial  profile  of  the  beam,  measured  1.5  m 
downstream  from  the  OKDE  interaction  region,  permits  a  measure  of 
the  ODKE  momentum  spread  We  will  discuss  each  of  the  various 
components  of  the  experiment  separately,  followed  by  our 
experimental  results. 

A  diagram  of  the  sodium  oven  is  shown  in  Figure  13.  The  oven 
consists  of  a  2.0  cm  diameter  cylinder  with  a  1-cm-long 
cylindrical  snout.  The  snout  has  an  inner  diameter  of  0.03  cm 
A  5.0-fim  wide  translatable  vertical  slit  is  located  12.5  cm  from 
the  center  of  the  oven  to  provide  initial  beam  steering  and 
collimation,  and  is  followed  by  a  variable-width  slit  127.5  cm 
from  the  oven  center,  and  a  fixed  25-fim  slit  just  before  the 
standing  wave  region  (165.5  cm  from  the  oven  center)  With 
careful  alignment  of  the  slits,  we  were  able  to  achieve  a  beam 
width  of  175  fim  at  the  spatial  beam  profile  monitor  located  at 
the  analysis  point  123  cm  from  the  standing  wave  interaction 
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Figure  12.  Optical  Kapi tza-Dirac  Effect:  experimental 
block  diagram. 
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Figure  13.  Sodium  oven  block  diagram. 
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Since  we  were  using  the  beam  spatial  profile  to  infer  the 
beam  momentum  distribution,  it  was  important  that  the 
longitudinal  velocities  of  the  atoms  be  uniform.  We  used  the 
technique  of  Steel  and  McFarlane,11  which  is  illustrated  in 
Figure  14,  to  prepare  the  beam  so  that  those  atoms  in  the  F=2 
ground  state  were  within  a  narrow  velocity  group.  Figure  15  is  a 
diagram  of  the  hyperfine  energy  levels  in  the  sodium  D2  line  at 
580  nm .  We  initially  pumped  the  atoms  with  light  that  was 
resonant  with  the  transition  from  the  F=2  ground  state  to  the  F=2 
excited  state  shown  in  the  figure  (beam  "A") .  The  laser  beam  was 
oriented  90*  to  the  beam,  ensuring  a  zero  Doppler  shift.  The 
atoms  in  the  F=2  excited  state  could  spontaneously  decay  to 
either  the  F=1  or  F=2  ground  states.  Since  atoms  in  the  F=1 
ground  state  would  no  longer  be  resonant  with  the  laser,  the 
number  of  atoms  in  the  F=2  state  would  eventually  be  negligible. 

The  second  step  of  the  velocity  selection  process  is  to  pass 
the  atoms  through  a  laser  beam  at  an  angle  of  approximately  45* 
to  the  atomic  beam  (laser  beam  "B"  in  the  figure).  If  the  atomic 
wave  vector  is  k  and  the  atomic  velocity  v,  the  laser  frequency 
is  Doppler  shifted  by  an  amount  -k*v  in  the  frame  of  each  atom. 
The  laser  frequency  was  chosen  to  interact  with  a  group  of  atoms 
having  a  particular  velocity  so  that  the  Doppler-shifted  light 
would  be  resonant  with  the  F=1  ->  F=2  transition  of  the  sodium  D2 
line.  These  excited  atoms  could  spontaneously  decay  to  either 
the  F=1  or  F=2  ground  state.  At  the  end  of  this  two-step 
process,  the  only  atoms  left  in  the  F=2  ground  state  were  those 
atoms  whose  velocities  were  within  the  narrow  band,  permitting 
them  to  be  resonant  with  the  45*  laser  beam.  We  were  able  to 
achieve  a  longitudinal  velocity  homogeneity  for  these  atoms  of 
Av/v=0.06.  The  rest  of  the  experiment  was  done  with  light  *^hat 
was  resonant  with  the  F=2  to  F=3  transition,  prohibiting  dipole 
decay  to  the  F=1  ground  state  and  hence  optical  pumping.  Atoms 
with  the  "wrong”  velocity  will  not  be  resonant  with  the  light, 
and  won’t  participate  in  subsequent  measurements. 


Figure  14.  Sodium  velocity  selection  scheme.  Left,  atomic 
ent  gy  level  scheme  showing  the  optical 
pumping  sequence  for  velocity  selection.  Right, 
geometry  for  the  optical  pumping  laser  beam  and 
the  velocity  selection  laser  beam. 
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The  next  stage  in  the  atomic  state  preparation  was  to 
experimentally  realize  a  strictly  two  level  system,  which  was 
done  by  pumping  the  atoms  with  circularly  polarized  light 
resonant  with  the  F=2  to  F=3  transition  A  small  magnetic  field 
of  1  to  3  gauss  was  applied  to  maintain  the  airection  of  the 
induced  atomic  dipole,  without  seriously  removing  the  degeneracy 
of  the  magnetic  substates.  The  selection  rule  for  excitation  is 
AmF=+l,  and  for  spontaneous  decay,  AmF=*l,0.  Thus  all  the  atoms 
traveling  into  the  interaction  region  are  left  in  the  mF  =  +  2 
ground  state.  Light  in  the  standing  wave  region  was  also 
circularly  polarized,  making  the  only  allowed  transition  one  from 
the  F=2,  mF  =  +  2  ground  state  to  the  F=3,  mF=+3  excited  state. 

This  process  is  illustrated  in  Figure  16 

The  OKDE  takes  place  in  the  standing  wave  region  shown  in 
Figure  17.  As  mentioned  in  the  preceding  paragraph,  the  standing 
wave  is  prepared  with  circularly  polarized  light  so  that  the 
strictly  two-level  characteristic  of  the  system  can  be  preserved. 
The  precise  standing  wave  alignment  is  verified  by  monitoring  the 
frequency  dependence  of  the  atomic  fluorescence  as  a  retro- 
reflecting  mirror  is  adjusted  to  reflect  the  initial  forward 
traveling  wave  back  onto  itself.  For  gross  misalignments  of  the 
two  beams,  as  the  laser  frequency  is  tuned,  a  double-peaked 
pattern  is  seen  in  the  fluorescence  spectrum  corresponding  to  the 
fluorescence  from  each  of  the  two  beams.  The  center  frequencies 
are  shifted  relative  to  each  other  because  of  the  slightly 
different  angles  (and  hence  Doppler  shifts)  that  the  two  optical 
beams  make  with  the  atomic  beam.  The  retro- ref  1 ected  beam  angle 
is  then  adjusted  so  that  the  two  fluorescence  peaks  "xactly 
overlap  each  other  in  frequency.  The  X/4  plate  and  polarizer, 
shown  in  Figure  17,  provide  for  circular  polarization  of  the 
light  and  for  optical  isolation  of  the  re tro- ref  1 ected  laser  beam 
from  the  dye  laser  With  careful  alignment,  we  were  able  to 
obtain  focused  beam  spot  sizes  of  7^m  FWHM .  This  small  spot  size 


F  i  gure 


Optical  pumping  mechanism  between  the  F  2  ground 
state  and  F=3  excited  state  in  the  D2  transition 
in  atomic  sodium  In  the  transition  shown  in  the 
figure,  atoms  i n  the  1  state  are  promoted  to 

the  mf -0  excited  state  by  a  photon  with  circular 
polarization  The  atom  can  then  decay  to  any  of 
the  ground  states  Eventually,  this 

mechanism  will  result  in  all  of  the  atoms  being 
in  the  mF-+2  ground  state 


Figure  17.  OKDE  interaction  region.  Retroref lected  light  is 
blocked  from  returning  to  the  laser  because  of  the 
90*  rotation  of  the  light  by  double  passing  the  X/4 
plate.  The  return  light  is  then  blocked  by  the 
linear  polarizer. 


ensures  that  atoms  will  traverse  the  OKDE  interaction  region  in 
less  than  one  spontaneous  decay  lifetime.  Spontaneous  decay 
would  have  introduced  a  diffusive  effect  on  the  spatial, 
distribution  which  would  have  obscured  the  diffractive  OKDE 
effect  we  were  studying. 

The  transverse  velocity  distribution  produced  by  the  OKDE 
results  in  a  spatial  distribution  if  the  atomic  beam  is  detected 
sufficiently  downstream  from  the  interaction  region.  The 
apparatus  we  developed  for  monitoring  the  spatial  profile  of  the 
beam  is  illustrated  in  Figure  18.  The  detector  is  essentially  a 
fluorescence  monitor  tuned  to  the  F=2  to  F=3  transition  so  that 
it  is  sensitive  only  to  the  atoms  within  the  velocity  group  of 
interest.  The  monitor  consists  of  a  laser  beam  that  is  reflected 
off  a  mirror  mounted  on  a  galvanometer.  The  galvanometer  is 
located  at  the  focus  of  a  cylindrical  lens  that  focuses  the  beam 
to  a  23  fim  spot.  The  accuracy  of  the  galvo-drive  circuit  is 
sufficient  to  reproducibly  position  the  spot  to  within  50  fim. 
Thus,  the  beam  effects  are  limited  more  by  the  175  fim  width  of 
the  atomic  beam  than  by  the  precision  of  the  spatial  profile 
monitor . 

Figure  19  is  the  spatial  profile  of  an  atomic  beam  before  and 
after  passing  through  the  standing  wave  region.  The  original 
beam  has  a  single  peak  that  is  split  into  two  peaks  that  are 
characteristic  of  the  theoretical  distributions  shown  in 
Figure  6.  The  experimental  resolution  was  insufficient  to 
resolve  the  recoil  due  to  single  photons,  which  would  have 
produced  the  comb  effect  seen  in  the  fully  quantum  mechanical 
calculation . 

Our  analysis  of  the  OKDE  data  has  been  enhanced  by  the 
development  of  two  computer  programs  that  run  on  the  HRL  VAX 
computers.  The  first  program  calculates  the  expected  spatial 
beam  profile  monitor  intensity  as  a  function  of  position, 
standing  wave  intensity,  standing  wave  spot  size,  and  atomic 
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Figure  19.  Experimental  data  showing  the  effect  of  the  internal 
standing  wave  on  the  spatial  profile  of  the  atomic 
beam.  The  noise  in  the  data  and  residual  bumps  is 
the  result  of  a  sodium  oven  malfunction  during  the 
experiment . 
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beam  size.  This  program  has  been  very  useful  in  understanding 
the  fundamental  limits  of  the  measurement  system,  in  discovering 
how  different  misalignments  affect  the  results,  and  in 
determining  what  kinds  of  signals  we  should  be  looking  for. 

The  second  program  is  used  for  analyzing  the  data  from  our 
TN-1710  signal  averager  on  the  VAX.  This  program  can  find  peaks 
and  widths  of  arbitrary  line  shapes,  calculate  statistical 
uncertainties  and  correlations  of  fit  parameters,  and  generate 
graphical  displays  of  the  data  and  fits  on  a  number  of  CRT  and 
hard  copy  output  units. 


SECTION  4 


ATOMIC  HYDROGEN  TWO-PHOTON  IONIZATION 


A.  INTRODUCTION 

One  of  the  factors  that  could  limit  the  maximum  cooling  rate 
of  an  atomic  hydrogen  beam  is  the  maximum  laser  beam  intensity 
possible  before  a  significant  number  of  hydrogen  atoms  are 
ionized  by  the  cooling  beam.  We  have  undertaken  experimental  and 
theoretical  studies  to  calculate  and  measure  the  magnitude  and 
the  features  of  this  effect. 

B.  THEORETICAL  STUDIES 

Multiphoton  ionization  of  hydrogenic  systems  is  concerned 
with  the  dynamics  involved  in  the  coupling  of  external  radiation 
fields  with  bound  quantum  systems  leading  to  the  production  of 
free  electrons  and  ions.  As  such,  multiphoton  ionization  has 
been  the  subject  of  intense  investigation  during  the  past  decade. 
One  of  the  simplest  and  most  fundamental  problems  concerns  the 
ionization  process  when  the  input  radiation  fields  induce  a 
resonant  two-photon  transition  into  the  continuum  from  the  ground 
state  of  a  hydrogenic  system.  In  this  case  there  are  two 
regimes.  If  the  radiation  is  not  resonant  with  a  bound-bound 
transition,  the  ionization  rate  can  be  computed  by  perturbation 
theory,  and  the  results  indicate  that  rate  is  proportional  to  the 
square  of  the  input  intensity.  However,  if  the  radiation  is 
resonant  with  a  bound-bound  transition,  saturation  processes  in 
the  population  difference  yield  an  expression  for  the  ionization 
rate  that  exhibits  saturation  behavior,  as  will  be  shown  below. 

Because  of  the  importance  of  using  Lyman -a  radiation  in  the 
cooling  and  trapping  of  atomic  hydrogen,  we  undertook  a  detailed 
description  of  the  role  of  ionization  of  hydrogen  atoms.  The 
objective  of  this  work  was  to  present  a  self-consistent 
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calculation  of  the  effects  of  resonance  radiation  on  the  two- 
photon  ionization  process  in  hydrogen  atoms.  The  results  showed 
that  two-photon  ionization  of  hydrogen  remains  negligible  for 
intensities  up  to  one  kw/cm2 . 

Above  10  kW/cm2,  significant  ionization  will  take  place.  The 
generation  of  such  a  plasma  environment  will  necessarily  increase 
the  emittance  of  the  beam,  possibly  leading  to  a  catastrophic 
breakup  as  it  propagates  in  space. 

1.  Formulation  of  the  Problem 

The  starting  point  in  our  analysis  was  the  formulation  of  a 
set  of  density  matrix  equations  that  describe  the  model  shown  in 
Figure  20.  The  bound  states  are  labeled  by  I i>  (i=l,2),  and  the 
continuum  by  the  set  {lk>}.  The  resonance  radiation  oscillates 
at  frequency  wL .  The  spontaneous  decay  rate  is  given  by  7.  In 
order  to  derive  a  set  of  equations  describing  the  evolution  of 
each  energy  level  via  the  density  matrix  formulation,  it  is 
important  to  understand  the  coherent  coupling  between  state  I 2> 
and  the  continuum  I k> .  This  understanding  comes  from  examining 
the  temporal  evolution  of  the  probability  amplitude  in  state  I 2> . 

Consider  the  Schroedinger  equation 

ih  ai*>/0t  =  (  H0  -  n  E  )  l*>  ,  (4.1) 

and  write  the  state  vector  !♦>  as  a  coherent  linear  superposition 
of  the  eigenstates  of  H0 , 

!♦>  =  cx(t)  ll>  exp(-iwLt)  +  c2(t)  I2>  exp(-iwLt) 

+  /  dk  c(k,t)  lk>  exp(-iwkt)  (4.2) 

Using  expression  (4.2)  in  Eq.  (4.1),  we  find  that  the  amplitudes 
c2  and  c(k,t)  satisfy  the  following  equations  of  motion; 
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Figure  20. 


Model  used  in  two-photon  ionization  calculations. 

The  two  bound  atomic  states  are  labeled  |1>  and  | 2> . 
The  continuum  states  are  labeled  |k>.  The  laser  is 
detuned  from  the  resonance  energy  E2-E!  by  an  amount 
A.  The  spontaneous  decay  rate  from  | 2>  to  |1>  is  7. 


ih  5  c2(t)/dt  =  -H21  E  cj(t)  exp(iu/21t) 

-/  dk  fi2k  E  c  (k ,  t)  exp(-iwk2t)  (4.3) 

ih  9  c(k,t)/9t  =  -fik2  E  c2(t)  exp(iwk2t)  ,  (4.4) 

where  pap  and  wap  are  the  electric  dipole  matrix  element  and 
frequency  separation  between  states  la>  and  !/?>,  respectively. 

We  assume  that  the  coupling  among  continuum  states  is  negligible 
and  that  the  off -resonant  two-photon  transition  between  the 
ground  state  and  the  continuum  has  a  small  contribution  in 
Comparison  with  the  resonant  transition. 

Equations  (4.3)  and  (4.4)  can  be  reduced  further  by  formally 
integrating  Eq.  (4.4)  and  using  the  result  in  Eq.  (4.3).  This 
procedure  leads  to  the  following  expression  for  c2 (t) : 

ih  9c2(t)/9t  =  -/i21  E  ct  (t)  exp(iw21t) 

+  (ifc)-1  /  dk  I  dt’  n2k  E(t)  Mk2  E(t’) 
x  exp[-iwk2(t-t')]  c2 (t  * )  .  (4.5) 

The  evolution  of  c2 (t)  in  the  presence  of  the  continuum  can 
be  deduced  in  a  simple  manner.  The  excitation  of  an  electron 
from  state  I 2>  by  the  radiation  field  can  only  take  place  if  t’ 
t.  This  simultaneity  requirement  implies  that  the  second  term  i 
Eq.  (4.5)  can  be  reduced  to  a  simple  decay  term  of  the  form  - 
rc2(t) .  T  =  ffl  is  the  intensity-dependent  ionization  rate  from 
the  excited  state  I 2> ,  and  a  is  the  single-photon  ionization 
cross  section. 

In  light  of  this  derivation  for  the  temporal  evolution  of 
c2(t),  the  density  matrix  equations  can  be  written  as 


a^n/at  =  'ip 2 2  +  i  R<>  P21  -  i  Ro*  P12 


(4.6) 


dp22/dt  =  p22  +  i  Ro*  P12  -  i  Ro  P 21  (4-7) 

ap12/at  =  i(w21  +  i7/2 )pl2  +  i  R0  (^22  -  Pn)  ,  (4.8) 

I 
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where  R0  is  the  on-resonance  Rabi  frequency.  The  set  of 
Eqs .  (4.6)  through  (4.8)  describes  the  evolution  of  the  bound 

system  in  the  presence  of  the  radiation  fields. 

2.  Exact  Analytical  Solutions 

Equations  (4.6)  through  (4.8)  have  an  exact  analytical 
solution  if  the  following  assumptions  are  made: 


(1)  The  radiation  field  has  a  pulse  duration  longer  than 
the  response  time  of  the  medium  (long  pulse 
approximation) . 

(2)  The  population  difference  between  the  bound  states  is 
slowly  varying  compared  to  the  smaller  of  the  natural 
linewidth  and  laser  detuning  from  intermediate  state 
resonance  (rate  equation  approximation) . 

(3)  The  laser  detuning  from  resonance  is  small  enough  so 
that  one  can  neglect  the  Bloch-Siegert  shifts  (rotating 
wave  approximation) . 


The  use  of  these  assumptions  in  the  formal  integration  of 
Eq.  (4.8)  gives  the  following  expression  for  the  optical 
coherence : 

P 12  =  i  R0  Ul2  (  P22  ~  P 11  )  >  (4-9) 

where  L12  =  {  iA  +  j/2  }  is  the  complex  Lorentzian  factor. 

Using  Eq.  (4.9)  in  Eqs.  (4.6)  and  (4.7),  we  obtain  a  closed 
system  of  equations  for  the  population  of  each  energy  level : 
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dpxl/dt  =  7p22  +  I R0  I  2  (p22  -  Pn)  Im(L12) 


(4.10) 


dp22/0t  =-(7+aI)p22  -  IR0I2  (p22  -  Pii)Im(L12)  ,  (4.11) 

where  Im(L12)  denotes  the  imaginary  part  of  L12.  We  shall  assume 
that  the  initial  conditions  of  the  problem  are  such  that 

Pi i (t=-«)  =  N0 

and 

P 22  (  t=-»  )  =  0. 

Nc  is  the  initial  density  of  the  atoms  in  the  ground  state. 

The  reduced  set  of  Eqs .  (4.10)  and  (4.11)  has  the  following 

exact  analytical  solutions: 

P72  -  {N0  K/X2>  exp  [-0 . 5 (7+aI+K) t]  sinh  (X  2 1)  (4.12) 

P  ii  =  Nj  -  p22  -  o"L  f  dt  p22(t  )  >  (4.13) 

where  K  =  I R0 I  2  Im(L12)  and  X2  =  {(7+aI+K)2  -  4aIK}1/2. 

Equations  (4.12)  and  (4.13)  constitute  the  main  result  of 
this  section  and  they  are  valid  in  the  regime  where  the  three 
fundamental  assumptions  (long  pulse,  rate  equation,  and  rotating 
wave  approximations)  are  valid. 

3.  Discussion 

Consider  hydrogen  atoms  excited  by  Lyman-a  radiation  tuned  to 
the  ls-2p  transition.  We  will  assume  a  cross  section  of  the 
order  of  10' 18  cm2  for  the  ionization  of  an  electron  located  in 
the  2p  state.  Figures  1  and  2  are  plots  of  the  temporal 
evolution  of  the  total  density  of  atoms  (i.e.,  pxl  +  p22). 


Figure  21  corresponds  to  the  case  where  the  Lyman-a  radiation  is 
tuned  on  resonance  to  the  ls-2p  transition,  and  shows  that 
negligible  ionization  takes  place  for  a  radiation  intensity  of 
1  kW/cm2  over  a  period  corresponding  to  100  times  the  spontaneous 
decay  lifetime  (1.6  ns).  As  the  intensity  is  increased  to 
500  kW/cm2,  signif  icant  ionization  takes  place,  leading  to  the 
possible  deterioration  of  the  beam  quality.  Figure  22  shows  the 
temporal  evolution  of  the  total  density  of  atoms  for  fixed 
intensity  as  a  function  of  the  radiation  detuning  from  line 
center .  Again,  for  zero  detuning,  significant  ionization  occurs. 
However,  as  the  detuning  is  set  to  10  7,  ionization  is  reduced 
significantly.  This  behavior  is  easily  explained.  If  the 
Lyman-a  is  on  resonance,  then  there  is  significant  excitation  of 
the  2p  state,  leading  to  a  stepwise  1 *sonantly  enhanced 
ionization  path.  However,  if  the  laser  is  tuned  off  resonance, 
negligible  excited  state  population  is  achieved  Hence,  the 
ionization  process  proceeds  along  a  non-resonant  pathway.  This 
explanation  is  illustrated  in  Figure  23,  where  the  excited  state 
population  is  plotted  as  a  function  of  time  As  is  shown  the 
excited  state  population  decreases  with  increasing  detuning  from 
resonance . 

A  quantity  of  interest  in  the  measurement  of  the  two-photon 
ionization  cross  section  is  the  two-photon  ionization  rate, 
defined  by, 


r  -  {9(Pii  *  P 2 2 ) /® t } / (p j j  *  p ; 2)  (4.14) 

Using  the  exact  analytical  expressions  for  p , j  and  p2?<  the 
resonant  two-photon  ionization  rate  is  given  by, 
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(aIK/X j ) /{coth (X j t)  ♦  (7+ctI+K) /2X 2 ) 


(4.15) 
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It  is  interesting  to  note  that  the  rate  is  time  dependent, 
because  the  exact  solutions  describe  the  temporal  behavior  of  the 
population  density  as  a  function  of  laser  intensity  and  detuning 
from  line  center.  The  long-time  behavior  (as  t  tends  to 
infinity)  is  given  by, 

r  =  2  0TK/(2X2+7+<7l+K)  ,  (4.16) 

which  is  a  nonlinear  function  of  the  laser  intensity.  Equation 
(4.16)  is  the  main  result  of  this  section;  it  reduces  to  the 
perturbation  solution  (proportional  to  I2)  in  the  low  intensity 
regime,  and  exhibits  saturation  behavior  when  the  single-photon 
ionization  rate  ffl  is  of  the  same  order  of  magnitude  as  the 
spontaneous  decay  rate  7 . 

4 .  Summary 

We  have  presented  a  theory  of  resonant  two-photon  ionization 
of  hydrogen  atoms.  We  have  obtained  exact  analytical  results  for 
the  total  density  of  the  atoms,  and  for  the  ionization  rate 
provided  under  the  three  assumptions  of  a  long  pulse,  the  rate 
equation,  and  the  rotating  wave  approximation.  Ionization  of  the 
hydrogen  atoms  can  be  avoided  if  the  Lyman-a  radiation  intensity 
is  restricted  to  less  than  1  kW/cm2 .  The  ionization  rate  was 
found  to  be  a  nonlinear  function  of  the  intensity  resulting  from 
saturation  of  the  bound  state  transition.  However,  the 
ionization  rate  reduces  to  the  result  obtained  from  perturbation 
theory  in  the  low  intensity  regime. 

C.  EXPERIMENTAL  STUDIES 

An  experimental  measurement  of  the  Lyman-a  two-photon 
ionization  cross  section  of  atomic  hydrogen  requires  development 
and  integration  of  both  a  tunable  narrow-band  coherent  photon 
source  and  a  high  current,  low  divergence  atomic  hydrogen  beam. 

We  h  ave  developed  the  hydrogen  beam  needed  for  the  experiment  and 
have  nearly  completed  the  optimization  of  the  Lyman-a  beam 
source.  Details  of  the  two  developments  all  given  below. 


1 .  Lyman-a  Production 


Figure  24  is  a  block  diagram  of  the  Lyman -a  source  we  have 
developed.  This  photon  source  consists  of  two  narrow  band 
oscillators  that  are  amplified  in  pulsed-dye  amplifiers  pumped  by 
the  outputs  of  a  dual-beam  Nd-YAG  laser.  The  outputs  from  the 
amplifiers  are  combined  in  a  mercury  cell  to  provide  the  desired 
Lyman -a  pulses. 

The  challenge  of  the  Lyman-a  system  operation  is  to  get 
concurrent  state-of-the-art  performance  from  the  many  complicated 
components  in  the  system.  Much  of  our  development  time  was  spent 
improving  the  reliability  of  the  coaaercikl  components  that  were 
used  in  the  system. 

The  narrow-band  oscillators  consist  of  two  Coherent  600-21 
ring  dye  lasers  pumped  by  two  Coherent  CR-18  argon  ion  lasers. 

The  first  oscillator  delivers  500  to  1000  mW  of  cw  light  at 
5454  A.  The  Rhodamine  560  dye  is  pumped  by  an  argon  ion  laser 
operating  at  6  W.  The  second  oscillator  contains  Kiton  Red  dye 
and  delivers  300  to  500  mW  of  power  at  6254  A.  The  Kiton  red  dye 
must  be  pumped  at  8  to  0  W. 

We  found  that  it  was  necessary  to  modify  the  Coherent  pump 
module  to  improve  the  lifetime  of  the  Rhodamine  560  dye  We  also 
developed  a  procedure  for  shutting  down  the  dye  lasers,  which 
resulted  in  an  order  of  magnitude  improvement  in  the  day -  to-  day 
reproducibility  of  our  dye  laser  outputs. 

The  outputs  of  the  dye  lasers  pass  through  optical  isolators 
consisting  of  a  pair  of  prisms  and  a  Faraday  rotator.  The  beams 
then  enter  modified  Quanta-Ray  PDA- 1  pulsed  dye  amplifiers  (PDAs) 
that  are  pumped  by  the  doubled  (532  nm)  outputs  of  a  dual -beam 
Nd-YAG  laser.  Required  modifications  in  the  PDAs  for  this 
application  included  replacing  or  modifying  of  the  internal 
telescoping  optics,  replacing  the  metal  pin  holes  with  diasK>nd 
ones,  and  mechanically  improving  the  pin  hole  mounts  After 
these  modifications,  we  were  able  to  get  10  ns  gaussian  pulses  in 
the  far  field  from  both  PDA’s. 
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Figure  24 


Experimental  setup  for  generation  of  VUV  Lyman-a 
rad  i at i on  The  output  from  the  dual  beam  Nd : YAG 

lasers  are  used  as  pumps  to  amplify  the  outputs 
of  the  tyo  narrow-band  oscillators  The  photons 
at  62$5  A  are  doubled,  yielding  photons  at 
3128  A  Light  at  3128  A  and  5455  A  is 
combined  in  an  He  cell  to  yield  a  sum  fre¬ 
quency  of  1216  A 


The  PDA  manufacturer  was  unable  to  supply  us  with  appropriate 
dye  concentrations  for  our  applications,  so  we  also  undertook  a 
study  to  learn  what  combinations  of  available  input  power,  pump 
power,  and  dye  concentrations  produced  the  best  results.  With 
optimum  concentrations  and  input  powers,  we  were  able  to  achieve 
20  mJ  at  5454  A,  and  6.5  mJ  at  6254  A.  The  radiation  at  6254  A 
was  doubled  in  a  KDP  crystal  to  produce  2.0  mJ  of  3127-A  light. 

The  radiation  at  3127  A  was  combined  with  the  5454-A  pulses 
in  the  mercury  cell  illustrated  in  Figure  25.  The  cell  consists 
of  a  center  region  that  is  heated  to  150*  to  200*C.  Mercury  from 
the  wick  evaporates  and  migrates  towards  the  end  regions  of  the 
cell,  which  are  water  cooled.  The  cooled  buffer  gas  in  the  end 
regions  condenses  the  mercury  back  onto  the  wicks.  The  mercury 
will  then  move  back  to  the  center  of  the  cell  where  the 
evaporation/condensation  cycle  begins  again.  Our  interaction 
region  is  approximately  10  cm  long.  We  verified  operation  of  the 
cell  by  measuring  the  absorption  of  light  from  the  mercury  lamp 
as  a  function  of  cell  temperature.  We  measured  an  aL  of  3.0  at  a 
cell  temperature  of  177*C.  After  several  weeks  of  operation,  we 
had  no  accumulation  of  mercury  on  the  cell  windows. 

The  mixing  process  in  the  mercury  cell  is  illustrated  in 
Figure  26.  Two  photons  at  3127  A  are  resonant  with  the  mercury 
transition  from  the  6s2  1S0  ground  state  to  the  6s7s  1S0  state. 
The  third  photon  at  5454  A  is  non-resonant  and  can  thus  be  tuned 
over  a  fairly  broad  range.  The  sum  frequency  generated  by  the 
addition  of  the  three  photons  is  at  1216  A  (Lyman-a) .  In  our 
tight  focusing  geometry,  the  maximum  signal  is  obtained  by  phase 
matching  the  sum  of  the  wave  vectors  of  the  input  and  output 
waves  so  that  bAk=-2 , 1 3  where  b  is  the  confocal  beam  parameter. 


We  have  used  two  diagnostics  to  study  the  radiation  exiting 
the  mercury  cell  —  a  nitric  oxide  (NO)  ionisation  chamber  and  a 
McPherson  Model  225  one  meter  scanning  VUV  monochromater . 

Figure  27  is  a  diagram  of  the  NO  ionization  detector.  The  cell 
was  filled  with  3.0  Torr  of  NO  gas,  which  ionizes  with  a  high 
efficiency  for  radiation  in  the  102-nm  to  135-nm  wavelength 
range.  With  our  geometry,  we  expected  an  efficiency  of  50% . 

With  the  NO  detector  attached  to  the  cell,  we  saw  a  small  ion 
signal  when  the  3127-A  light  was  introduced  to  the  mercury  cell. 
The  signal  was  resonant  with  the  two-photon  resonance  we  were 
tuning  for,  but  didn’t  show  a  strong  dependence  on  the  intensity 
of  the  5454-A  light.  Also,  the  tuning  width  was  broader  than  we 
would  have  expected  from  either  the  natural  line  width  of  the 
mercury  resonance  or  the  Doppler  broadening  of  the  atoms  in  the 
mixing  cell. 

To  better  understand  our  results,  we  attached  the  mercury 
cell  to  the  McPherson  monochromater.  The  monochromater  was  first 
calibrated  with  a  hydrogen  Hinterreger  light  source  at  Lyman-a 
and  with  a  number  of  lines  from  a  mercury  lamp.  The  radiation 
was  detected  with  a  solar  blind  photomultiplier  tube.  We  have 
not  yet  been  able  to  reproduce  on  the  monochromater  the  VUV 
signals  we  saw  in  the  NO  detector. 

2.  Hydrogen  Beam  Production 

Figure  28  is  a  diagram  of  our  hydrogen  beam  production 
facility.  Molecular  hydrogen  is  dissociated  by  an  rf  field  in  a 
water  cooled  glass  cell.  A  snout  at  the  end  of  the  cell  is 
capped  with  an  ef fuser  plate  consisting  of  a  close  packed  array 
of  2 . 0-/im-diameter ,  3.0-mm-long  capillaries.  These  capillaries 
provide  a  differential  pressure  between  the  rf  discharge  region 
and  the  high  vacuum  region,  and  they  provide  good  beam 
collimation.  With  some  experimentation,  we  were  able  to  achieve 
a  beam  divergence  of  less  than  1*.  At  this  divergence,  the 
natural  line  width  of  the  hydrogen  Lyman-a  transition  will  be 
comparable  to  the  Doppler  line  width. 


The  hydrogen  beam  intersects  the  Lyman-a  beam  in  the  first 
chamber.  Before  intersecting  the  hydrogen  beam,  howeverr,  the 
photon  beam  is  first  dispersed  by  a  diffraction  grating  to 
separate  out  photons  (which  also  exit  the  mercury  cell)  with 
unwanted  wavelengths.  A  hydrogen  discharge  cell  is  provided  in 
the  photon  beam  line  to  monitor  the  frequency  of  the  radiation. 
Radiation  flux  is  measured  with  an  NO  ionization  detector. 

The  hydrogen  flux  can  be  monitored  by  two  devices.  The  best 
absolute  measurement  is  given  by  a  bolometer,  consisting  of  a 
platinum  strip  that  intersects  the  entire  hydrogen  beam.  Atomic 
hydrogen  incident  on  the  strip  will  recombine  on  the  surface  of 
the  platinum  and  leave  the  surface  as  molecular  hydrogen.  The 
total  energies  of  recombining  hydrogen  atoms  will  heat  the  strip, 
thus  changing  its  resistance.  After  balancing  a  bridge  circuit 
that  has  the  platinum  strip  as  one  of  its  legs,  the  hydrogen  beam 
is  turned  off  and  the  platinum  strip  is  electrically  heated  until 
its  resistance  matches  that  of  the  strip  as  it  was  being  heated 
by  the  hydrogen  beam.  The  product  of  the  electrical  current  and 
voltage  supplied  to  the  strip  gives  the  power  that  must  also  have 
been  supplied  by  the  recombining  hydrogen  atoms  in  the  beam. 
Dividing  this  power  by  the  energy  of  recombination  of  the 
hydrogen  atoms  yields  the  number  of  hydrogen  molecules  leaving 
the  surface  of  the  platinum,  per  second. 

We  also  used  a  UTI  mass  spectrometer  as  a  secondary  monitor 
of  the  beam  flux.  Once  calibrated  against  the  platinum  strip 
detector,  this  spectrometer  provides  a  real-time  measure  of  the 
hydrogen  flux,  which  can  be  monitored  continuously  during  data 
collection.  We  were  able  to  generate  atomic  hydrogen  beam 
currents  of  6x10* 7  atoms/s. 


SECTION  6 


w -PULSE  COOLING 

A.  INTOODUCTI ON 

Experimental  studies  of  laser  modification  of  atomic 
velocities  using  the  phenomenon  of  laser- induced  fluorescence  as 
proposed  by  Hansch  and  Schawlow,14  have  been  quite  extensive. 
Evidence  of  the  slowing  of  an  atomic  beam  by  resonance  radiation 
pressure  was  demonstrated  by  Bjorkholm  et  al . , 1 6  Balykin  et 
al  .  ,  1  *  Phillips  and  Metcalf,17  and  Prodan  et  al . 1  *  Recently, 
experiments  by  Prodan  et  al . 1 9  and  Ertmer  et  al.2*  demonstrated 
that  an  effusive  atomic  beam  was  stopped  by  laser  light.  These 
works  were  followed  by  the  recent  experimental  demonstration  of 
radiation  trapping  of  atoms  by  Chu  et  al . 2 1 

In  the  Hansch  and  Schawlow  cooling  procedure,  the  atoms  are 
illuminated  with  counterpropagating  laser  beams  that  are  tuned 
into  the  upper  half  of  the  Doppler  contour  of  a  resonance  line, 
the  laser  frequency  either  covering  the  full  upper  half  of  the 
contour,  or  scanned  toward  line  center  as  cooling  progresses. 
Cooling  results  from  the  net  momentum  transfer  from  the  laser 
photons  to  the  atoms  during  resonance  fluorescence  events.  The 
cooling  rate  is  thus  limited  by  the  natural  decay  rate  of  the 
atoms  from  the  spontaneous  radiation;  the  maximu'  ooling  rate 
occurs  at  laser  power  levels  that  saturate  the  transition. 

The  above  limitation  on  the  cooling  rate  could  be  lifted  if 
it  were  possible  to  alternate  the  direction  and  detuning  of  the 
separate  upward-stimulating  and  downward-stimulating  photons. 
Then,  a  net  momentum  transfer  to  the  atom  could  be  made  to  occur 
at  transition  rates  greater  than  the  spontaneous  decay  rate. 

This  goal  can  be  achieved  by  using  a  train  of  alternately 
oppositely  directed  and  oppositely  detuned  n  pulses. 


The  concept  of  modifying  the  velocity  of  the  atoms  using  x 
pulses  was  advanced  by  Kazantsev  in  1974. 22  A  similar  scheme  was 
proposed  by  Nebenzahl  and  Szoke23  in  their  studies  of  atomic-beam 
deflection  using  stimulated  emission  processes.  Recently,  Aspect 
et  al . 3  demonstrated  the  collimation  of  a  cesium  atomic  beam 
using  stimulated  emission.  A  related  application  using  x  pulses 
to  deflect  an  atomic  beam  for  isotope  separation  was  considered 
by  Friedman  and  Wilson.24 

B.  PHYSICAL  PICTURE 

The  basic  physics  of  X  pulse  cooling  is  described  as  follows: 
consider  the  interaction  of  a  two-level  atom  with  a  pair  of  x 
pulses  incident  from  opposite  directions.  First,  a  photon  from 
the  x  pulse  from  the  right  is  absorbed,  exciting  the  atom  to  the 
upper  state.  An  increment  of  momentum  -hk  has  been  transferred 
to  the  atom.  Here  k  is  the  wavenumber  and  K  is  Planck’s  constant 
divided  by  2 x.  Next,  the  x  pulse  from  the  left  stimulates  the 
atom  to  emit  a  photon  to  the  left,  resulting  in  a  recoil  momentum 
transfer  to  the  atom  of  -hk .  Hence  the  net  momentum  transferred 
to  the  atom  from  the  interaction  of  both  x  pulses  is  -2hk . 
Clearly,  a  train  of  n  alternately  counter-propagating  x  pulses 
sequentially  incident  upon  the  atom  will  result  in  a  momentum 
transfer  of  nhk  in  the  direction  of  propagation  of  the  upward 
stimulating  x  pulses,  provided  that  the  duration  of  the  pulse 
train  remains  short  in  comparison  with  the  radiative  lifetime  of 
the  transition.  To  retard  the  motion  of  the  atom’s  the  upward 
stimulating  x  pulses  should  be  directed  opposite  to  the  motions 
of  the  atom  and  tuned  above  resonance,  while  the  downward 
stimulating  x  pulses  should  be  oppositely  directed  and  tuned 
below  resonance. 

To  cool  a  one-dimensional  velocity  distribution  of  atoms,  a 
train  of  x  pulses  with  a  detuning  sequence  chosen  to  retard  the 
motion  of  atoms  moving  in  one  direction  must  be  followed  by  a 
train  with  the  start  of  the  detuning  sequence  shifted  to  the 


oppositely  directed  beam  in  order  to  retard  the  motion  of  the 
atoms  moving  in  the  opposite  direction.  To  ensure  cooling  of  the 
distribution,  each  ir-pulse  train  must  effect  a  larger  momentum 
transfer  to  atoms  in  the  half  of  the  distribution  for  which  the 
motion  is  being  retarded,  than  to  the  other  half.  This  desired 
result  will  occur  if  the  frequency  of  the  ir-pulse  train  is  chosen 
to  be  resonant  with  atoms  with  Doppler  shift  of  half  the  Doppler 
width  of  the  distribution,  since  the  excitation  state  of  the 
atoms  on  that  side  of  the  distribution  will  remain  in  coherence 
with  the  w-pulse  train  for  a  larger  number  of  w  pulses  than  will 
atoms  in  the  other  side  of  the  distribution.  The  two  ir-pulse 
trains  must  be  separated  by  a  few  radiative  lifetimes  in  order  to 
reset  the  excitation  state  of  all  the  atoms  to  the  ground  state 
before  switching  the  detuning  sequence. 

This  physical  description  of  the  cooling  process  achieved 
with  the  ir-pulse  technique  is  rigorously  justified  using  quantum 
mechanical  transport  equations.  The  formulation  of  the  equations 
is  the  core  of  the  next  section. 

C.  THE  QUANTUM  TRANSPORT  EQUATIONS,  INCLUDING  ATOMIC  RECOIL 

The  quantum  transport  equations  were  first  derived  by 
Kolchenko  et  al?6  when  they  studied  the  spectral  response  of 
gaseous  systems  in  saturated  absorption  spectroscopy.  They 
observed  that  the  origin  of  the  symmetrical  splitting  of  the 
Doppler-free  absorption  line  was  the  atomic  recoil  experienced  by 
the  atom  upon  absorbing  and  emitting  photons  The  same  equations 
have  been  rediscovered  by  numerous  authors  This  section 
presents  a  detailed  derivation  of  the  quantum  transport 
equations,  which  take  into  account  the  effects  of  Doppler 
broadening,  grating  washout,  and  atomic  recoil 

If  R  denotes  the  center -of -mass  displacement  vector  and  r  the 
internal  coordinates,  then  the  space-time  evolution  of  an  atom  is 
governed  by  the  Schrodinger  equation, 


ifc  9  *(R,r,t)/9t  =  [-R2V2/2M  +  H0  +  V(R , t) ] * (R , r , t) , 


(5.1) 


where  the  first  term  on  the  right  hand  side  accounts  for  the 
translational  motion  of  the  atom,  H0  is  the  electronic 
Hamiltonian,  and  V(R,t)  is  the  interaction  Hamiltonian.  We  will 
consider  an  interaction  potential  for  the  induced  electric  dipole 
moment  of  an  atom  interacting  with  the  radiation  field. 

We  begin  by  expanding  the  wavefunction  f  as  a  linear 
superposition  of  eigenstates  <pm  of  H0  : 

*(R,r,t)  =  E  cw(R,t  )<pm  exp  ,  (5.2) 

where  the  probability  amplitude  cm  accounts  for  the  translational 
motion  of  the  atom  in  specific  eigenstate  .  Using  expression 
(5.2)  in  the  Schrodinger  equation  (5.1),  one  finds  the  temporal 
evolution  of  cn  is  described  by 

ih  9cn/9t  =  -h2V2  cn/2M  +  E  <nlV/m>  exp(iwmnt)  cm  (5.3) 

where  <nlVlm>  and  wnw  are  the  matrix  element  of  V  and  frequency 
difference  between  eigenstates  I m>  and  lm>,  respectively. 

Our  goal  is  to  obtain  an  evolution  equation  for  the 
observables  of  the  problem  in  question.  In  this  case,  the 
observables  are  the  population  of  each  eigenstate  and  the  optical 
coherence  generated  by  the  radiation  field.  Defining  the  Wigner 
distribution  function  as 


f„.(R.P.t)  =  (l/2*h) 3  /  dR’  exp (-ip .R ’ /R) 

Pnmi  R+R’/2,R-R72)  ,  (5.4) 

where  pnm  is  the  observable,  and  using  the  Wigner  distribution  in 
expression  (5.4),  we  find  that  the  quantum  mechanical  transport 
equation  is  described  by, 
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ih(a/dt  +  P*v/M)  fnai(R,p,t) 

=  E/dk  [<niv’la>  exp(-ik*R) fam  (R,p+Rk/2, t)exp(iw„at) 

-  f  „a(R,p-hk/2,t)<al  V' lm>exp(-ik»R)exp(iwai,t)]  ,  (5.5) 

where  the  convective  derivative  p»V/M  accounts  for  the 
translational  motion  of  the  atom  and  gives  rise  to  Doppler 
broadening  and  grating  washout.  The  quantity  <nlV’la>  is  the 
spatial  Fourier  transform  of  <nlVla>.  Equation  (5.5)  is  the  main 
result  of  this  section:  it  describes  exactly  the  interaction  of 
an  atom  with  a  classical  radiation  field,  and  takes  into  account 
the  atomic  recoil  effect. 

D.  NUMERICAL  SOLUTIONS 

In  essence,  Eq.  (5.5)  is  the  quantum  generalization  of  the 
classical  Boltzmann  equations.  Applying  it  to  the  case  of  two 
level  atoms,  we  can  show  that  Eq.  (5.5)  reduces  to  equations 
describing  the  following  distribution  functions:  M(p,t)  = 
f  2  2  +  fii>  the  number  density  of  atoms;  N(p,t)  =  f22  ~  fii»  the 
population  difference  between  the  two  states;  and  f  *  2  (p  >  t)  ,  the 
optical  coherence.  The  quantum  transport  equations  are 


a  N/at  =  4  n  f, (p,t) 

3  M/8t  =  2  fl  [Fj(p+Rk,t)  -  F j  (p , t) ] 
(d/at  +  iAw  -  ip*k)  F(p,t) 


(5.6) 


(5.7) 


=  i  (1  N(p,t)  +  i  fl  [M(p+hk,t)  -  M(p,t)],  (5.8) 

where  (l=pE/h  is  the  on-resonance  Rabi  frequency,  fl  is  the  dipole 
moment  of  the  transition,  and  E  is  the  electric  field  amplitude. 
We  have  made  the  substitution:  F  =  f12  exp  (ikx-iAwt) ,  where  Aw 
is  the  detuning  from  line  center.  F  •,  indicates  the  imaginary 
part  of  F. 
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We  have  carried  out  numerical  solution  of  Eqs  (56)  through 
(5.8),  using  a  train  of  w  pulses  as  the  applied  field,  with 
alternate  pulses  oppositely  directed  and  tuned  on  opposite  sides 
of  line  center,  as  required  for  cooling.  Hydrogen  atoms  were 
selected  for  the  atomic  system  in  question,  and  the  wavelength  of 
the  radiation  field  was  assumed  to  be  tuned  to  the  Lyman-a 
transition  of  hydrogen.  The  results  of  the  simulation  are 
detailed  in  Figure  29,  which  shows  an  initially  Maxwellian 
momentum  distribution,  and  the  distribution  after  the  application 
of  16  ir-pulses  tuned  in  the  manner  discussed  above  Also  shown 
are  the  plots  of  I  FI  after  the  application  of  each  of  the  » - 
pulse  trains  tuned  to  the  opposite  side  of  the  distribution  We 
can  see  from  these  plots  that  coherence  is  maintained  better  on 
the  side  of  the  distribution  that  is  being  cooled  than  't  is  on 
the  opposite  side,  as  expected  As  required,  the  phases  of  all 
of  the  atoms  were  reset  to  zero  before  the  application  of  the 
second  train  of  w-pulses  Also  evident  is  an  essentially 
symmetric  distribution,  M(p),  that  is  cooled,  except  at  the  wings 
of  the  distribution  where  some  apparent  heating  of  the 
distribution  has  occurred. 

B.  SUMMARY 

In  conclusion,  we  have  examined  radiation  cooling  of  a 
distribution  of  two-level  atoms  in  the  coherent  pulse  regime  A 
train  of  oppositely  directed  and  oppositely  detuned  e-pulses, 
with  specific  zero-field  periods  between  the  ir-pulses,  is 
required  for  optimal  cooling.  The  cooling  was  desKinstrated 
through  numerical  analysis  of  the  quantum  transport  equations 
The  scheme  proposed  here  provides  an  alternate  method  for  cooling 
neutral  particle  beams. 
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SECTION  6 


SUMMARY 

In  this  report,  we  have  summarized  the  work  completed  in 
three  areas  relevant  to  cooling  and  trapping  of  neutral  atomic 
beams.  We  first  provided  a  summary  of  ways  in  which  radiation 
forces  can  be  used  to  reduce  the  emittance  or  trap  neutral  atoms 
This  summary  was  followed  by  a  brief  review  of  the  physics  of 
atom/laser  interactions.  Finally,  we  reviewed  our  theoretical 
and  experimental  investigations  in  Optical  Kapi tza-Dirac  Effect, 
hydrogen  two- photon  resonant  ionization  at  Lyman-a,  and  ir  pulse 
cool ing . 

In  our  Optical  Kapi tza-Di rac  work,  we  observed  the 
diffraction  effect  predicted  by  theory.  As  a  practical  matter, 
the  effect  will  not  become  significant  in  the  time  and  intensity 
regimes  expected  for  cooling  relativistic  atomic  beams.  Our 
theoretical  studies  of  the  effects  of  very  intense  beams  on  the 
OKDE  were  also  presented.  Our  dressed  state  formalism  provides 
means  of  treating  the  effects  of  intense  fields  for  large 
detunings  and  intensities  above  the  saturation  intensity. 

The  atomic  two-photon  ionization  work  also  involved  both 
experimental  and  theoretical  investigations.  Our  calculations 
show  that  ionization  will  not  be  a  problem  in  cooling  a  hydrogen 
beam,  provided  the  intensity  is  kept  below  1  kW/cm2 .  Since  this 
figure  is  considerably  larger  than  the  Lyman-a  saturation 
intensity,  it  is  unlikely  this  limit  will  prove  to  be  a  serious 
one.  The  experimental  efforts  have  resulted  in  development  of  a 
hydrogen  beam  source  and  a  VUV  source  that  should  prove  adequate 
to  perform  the  indicated  measurements.  However,  further  work 
must  be  done  on  the  VUV  source  before  the  measurements  can  be 
completed . 

The  theory  of  K  pulse  cooling  was  developed  in  Section  5. 
This  technique  appears  to  be  a  viable  alternative  to  the  more 


traditional  approach  of  cooling  by  spontaneous  emission.  A 
numerical  simulation  of  cooling  established  the  feasibility  of 
this  alternative . 

In  conclusion,  laser  cooling  appears  to  be  an  attractive 
technique  for  cooling  and  trapping  neutral  atoms.  The  two 
limitations  imposed  by  standing  wave  diffraction  effects  and 
photo-ionization  are  interesting  and  real,  but  aren’t  expected  to 
seriously  affect  the  performance  of  present  cooling  and  trapping 
schemes.  Alternative  approaches  to  cooling  appear  to  offer 
advantages  over  traditional  cooling  techniques  that  rely  on 
spontaneous  emission. 
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The  concept  of  radiation  cooling  with  r  pulses  is  presented  and  quantified  with  a  theoretical  analysis.  Numerical 
analysis  of  e-pulse  cooling  on  the  hydrogen  Lyman-o  transition  is  presented  as  an  example.  The  cooling  rate  and 
limiting  temperature  are  discussed  and  compared  with  steady-state  radiation  cooling. 


1.  INTRODUCTION 

The  possibility  of  using  light  to  modify  the  velocity  distribu¬ 
tion  of  an  ensemble  of  atoms  appears  to  have  originated  with 
the  pioneering  work  of  Kolchenko  et  a/.,1  who  analyzed  the 
effect  of  atomic  recoil  in  the  spectral  line  shape  of  an  atom 
interacting  with  a  resonance  radiation  field.  At  almost  the 
same  time,  Ashkin2  proposed  the  concept  of  resonantly  en¬ 
hanced  radiation  pressure  in  a  gaseous  system.  These  works 
stimulated  the  development  of  two  important  contributions 
in  the  area  of  laser  cooling.  First  was  the  concept  of  laser- 
induced  cooling  of  gases  by  Hansch  and  Sc  haw  low.1  which  is 
the  foundation  of  much  of  today's  research  in  the  area. 
Second  was  the  concept  of  atomic-beam  birefringence  under 
the  action  of  a  strong  radiation  field  by  Kazantsev.4  which 
was  the  cornerstone  of  the  so-called  optical  Stern-Gerlach 
effect.  The  latter  was  recently  demonstrated  by  Moskowitz 
et  a/..'  and  the  theoretical  analysis  of  Kazantsev  was  con¬ 
firmed  by  Cook  and  Bernhardt.6 

Experimental  studies  of  laser  modification  of  atomic  ve¬ 
locities  that  use  the  phenomenon  of  fluorescence  have  been 
quite  extensive.  Evidence  of  the  slowing  of  an  atomic  beam 
by  resonance  radiation  pressure  was  demonstrated  by  Bjork  - 
holm  et  al.'  Balykin  et  a/.,8  Phillips  and  Metcalf,8  and  Pro- 
dan  et  al.10  Recently,  experiments  by  Prodan  ef  a/.11  and 
Ertmer  et  al.12  demonstrated  that  an  effusive  atomic  beam 
was  stopped  by  laser  light.  These  works  were  followed  by 
the  recent  experimental  demonstration  of  radiation  trap¬ 
ping  of  atoms  by  Chu  et  al 13 

In  the  Hiinsch-Schawlow  radiation  cooling  procedure,  the 
atoms  are  illuminated  with  counterpropagating  laser  beams 
that  are  tuned  into  the  upper  half  of  the  Doppler  contour  of  a 
resonance  line,  the  laser  frequency  either  covering  the  full 
upper  half  of  the  contour  or  scanned  toward  line  center  as 
cooling  progresses.  Cooling  results  from  the  net  momentum 
transfer  to  the  atoms  from  the  laser  photons  during  reso¬ 
nance  fluorescence  events.  The  cooling  rate  is  thus  limited 
by  the  natural  decay  rate  of  the  atom  due  to  spontaneous 
radiation,  the  maximum  cooling  rate  occurring  at  laser  pow¬ 
er  levels  that  saturate  the  transition.  An  increase  of  laser 
power  above  saturation  does  not  increase  the  cooling  rate 
since  absorption  followed  by  stimulated  emission  from  the 
same  beams  does  not  transfer  any  net  momentum  to  the 
atom. 

The  above  limitation  on  the  cooling  rate  could  be  lifted  if 
it  were  possible  to  alternate  the  direction  and  detuning  of 


the  separate  upward-stimulating  and  downward-stimulat¬ 
ing  photons.  Then,  a  net  momentum  transfer  to  the  atom 
could  be  made  to  occur  at  transition  rates  greater  than  the 
spontaneous  decay  rate.  As  will  be  illustrated  below,  this 
can  be  accomplished  with  the  use  of  a  train  of  alternately 
oppositely  directed  and  oppositely  detuned  »  pulses.  The 
concept  of  modifying  the  velocity  of  atoms  using  r  pulses 
was  advanced  by  Kazantsev14  in  1974.  A  similar  scheme  was 
proposed  by  Nebenzahl  and  Szoke15  in  their  studies  of 
atomic-beam  deflection  using  stimulated-emission  process¬ 
es.  A  related  application  of  using  *  pulses  to  deflect  an 
atomic  beam  for  the  purpose  of  isotope  separation  was  con¬ 
sidered  by  Friedman  and  Wilson.16 

The  objective  of  this  paper  is  to  identify  the  critical  issues 
associated  with  ir-pulse  cooling  of  a  distribution  of  atomic 
velocities,  to  estimate  the  cooling  rate  and  limiting  tempera¬ 
ture.  and  to  compare  these  performance  figures  with  those 
for  steady-state  laser  cooling.17  Below.  Section  2  describes 
the  physical  picture  of  e-pulse  cooling  and  the  restrictions 
on  the  incoming  tram  of  *  pulses  due  to  Doppler  dephasing. 
A  rigorous  theoretical  foundation  for  ir-pulse  cooling  that 
utilizes  the  Wigner  distribution  function  is  presented  in 
Section  3.  The  results  of  a  numerical  integration  of  the 
quantum-mechanical  transport  equations  applied  to  ir- 
pulse  cooling  on  the  hydrogen  Lyman-a  transition  appear  in 
Section  4.  Section  5  compares  the  cooling  rates  and  re¬ 
quired  power  levels  for  steady-state  and  »- pulse  cooling,  and 
Section  6  summarizes  our  results. 

2.  PHYSICAL  PICTURE 

The  basic  physics  of  ir-pulse  cooling  is  as  follows:  Consider 
the  interaction  of  a  two-level  atom  with  a  pair  of  ir  pulses 
incident  from  opposite  directions.  First,  a  photon  from  the 
ir  pulse  from  the  right  is  absorbed,  exciting  the  atom  to  the 
upper  state.  An  increment  of  momentum  -Ihk)  has  been 
transferred  to  the  atom,  i  Here  k  is  the  photon  wave  number 
and  h  is  Planck's  constant  divided  by  2».)  Next,  the  <r  pulse 
from  the  left  stimulates  the  atom  to  emit  a  photon  to  the  left, 
resulting  in  a  recoil  momentum  transfer  to  the  atom  of 
-(hfcl.  The  net  momentum  transferred  to  the  atom  due  to 
the  interaction  with  both  ir  pulses  is  —  <  2A A I .  Clearly,  a  train 
of  ,V  alternately  counterpropagating  ir  pulses  sequentially 
incident  upon  the  atom  will  result  in  a  momentum  transfer 
of  Sthk)  to  the  atom  in  the  direction  of  propagation  of  the 
upward-stimulating  <r  pulses,  provided  that  the  duration  of 
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the  pulse  train  remains  short  in  comparison  with  the  radia 
live  lifetime  of  the  transition.  To  retard  the  atom's  motion, 
the  upward -stimulating  r  pulses  should  thus  be  directed 
opposite  to  the  atom's  motion  and  tuned  above  resonance, 
while  the  downward-stimulating  pulses  should  be  oppositely 
directed  and  tuned  below  resonance 

To  cool  a  one-dimensional  velocity  distribution  of  atoms, 
a  train  of  »  pulses  with  a  detuning  sequence  chosen  to  retard 
the  motion  of  atoms  moving  in  one  direction  must  be  fol¬ 
lowed  by  a  train  with  the  start  of  the  detuning  sequence 
shifted  to  the  oppositely  directed  beam  in  order  to  retard  the 
motion  of  the  atoms  moving  in  the  opposite  direction  To 
ensure  cooling  of  the  distribution,  each  x-pulse  train  must, 
of  course,  effect  a  larger  momentum  transfer  to  atoms  in  the 
half  of  the  distribution  for  which  the  motion  is  being  retard¬ 
ed  than  to  the  other  half  This  will  occur  if  the  frequency  >f 
the  e-pulse  train  is  chosen  to  be  resonant  with  atoms  with  a 
Doppler  shift  of.  say.  half  the  Doppler  width  of  the  distribu¬ 
tion.  since  the  excitation  state  of  the  atoms  in  that  side  of  the 
distribution  will  remain  in  coherence  with  the  r -pulse  train 
for  a  larger  number  of  *  pulses  than  will  atoms  in  the  other 
side  of  the  distribution.  The  two  *  pulse  trains  must  be 
separated  by  a  few  radiative  lifetimes  in  order  to  reset  the 
excitation  state  of  all  the  atoms  to  the  ground  state  before 
switching  the  detuning  sequence 

It  is  interesting  to  inquire  where  the  entropy  of  the  atom 
motion  goes  in  the  case  of  cooling  with  a  train  ot  ir  pulses  In 
the  case  of  steady-state  laser  cooling,  the  entropy  is  carried 
away  by  the  spontaneously  emitted  photons  During  cool¬ 
ing  with  a  train  of  r  pulses,  spontaneous  emission  is  assumed 
not  to  occur.  However,  the  entropy  generation  still  appears 
in  the  radiated  photons.  Radiation  from  the  atoms  during 
illumination  by  a  e-pulse  train  is  accounted  for  by  the  in¬ 
duced  transition  dipole  moment,  which  is  proportional  to 
the  off-diagonal  terms  in  the  density  matrix.  F in  isee  be 
low  i  The  generated  entropy  comes  from  nonresonant 
atoms'  emitting  photons  into  a  larger  number  of  temporal 
modes  than  do  the  resonant  atoms.  As  discussed  above,  the 
more  randomized  phase  of  Fin  for  the  nonresonant 
atoms  also  accounts  for  the  loss  of  phase  sequencing  with  the 
» -pulse  train  necessary  for  adjacent  *  pulses  to  produce 
equally  directed  forces  on  the  atom  Since  Fin  is  cycling 
much  faster  than  the  spontaneous  emission  rate,  the  rate  of 
entropy  generation  can  be  larger  than  in  the  case  of  steady- 
state  laser  cooling,  thus  permitting  a  faster  cooling  rate  for 
r  pulse  cooling 

The  accumulated  rate  of  phase  error  between  the  excita 
tion  state  and  the  r-pulae  train  that  occurs  for  nonresonant 
atoms  will  determine  the  maximum  number  of  r  pulses  that 
can  be  used  in  a  given  e-pulse  train  Also,  the  greater  the 
disparity  of  the  accumulated  phase  error  on  the  different 
sides  of  the  distribution,  the  faster  the  cooling  rate  w  ill  be  for 
a  given  w  pulse  It  turns  out  that  the  accumulated  phase 
error  can  be  minimized  and  the  disparity  of  phase  error 
maximized  through  the  use  of  s  specific  zero  field  period 
between  the  *  pulses  This  technique  is  best  appreciated 
through  reference  to  the  density-matrix  vector  model  for  the 
evolution  of  the  coherent  excitation  state  of  a  two  level 
atom.1* IS 

Figure  l  shoves  the  density  matrix  vector  for  one  atomic 
velocity  before  the  application  of  the  first  *  pulse,  after  the 
zero-field  period,  and  after  the  application  of  the  second  » 


Fig  I  Vector  representation  of  the  state  ol  a  two  level  atom  lor 
one  atomii  velocity  for  time  zero  alter  the  first  »  pulse  r  after  the 
zero  tield  perusl  '  and  alter  the  setond  r  pulse  : , 

pulse  The  density  matrix  vector  rotates  about  the  vector  B 
“  il.’ei  +  Awce^i.  where  1?  »  UE  h  is  the  Rabi  frequency  u  is 
the  transition  dipole  moment.  E  is  the  laser  field,  and  A*.  is 
the  detuning  from  resonance  It  rotates  about  B  at  an  angu 
lar  rate  equal  to  1 11- +  A-.1  -  When  the  vector  points  in  the 
negative  e >  direction  the  atom  is  in  the  ground  state  w  hen  it 
points  in  the  positive  e  direction  the  atom  is  in  the  upper 
state  All  other  directions  represent  the  atom  in  a  superpo 
sit  ion  of  these  two  states 

It  can  be  seen  from  this  picture  that  if  there  is  no  field  off 
period  between  the  *  pulses,  then  the  phase  error  spread 
accumulates  to  within  '-<1.’  A-.  i  pulses  However  dur 
mg  a  zero  field  interval  between  r  pulses  the  vector  rotates 
about  the  e.  axis  at  a  rate  equal  to  the  detuning  Thus,  it  we 
choose  the  duration  ot  the  off  field  interval  to  be  r  <  A- ;  o 
where  A-/,  is  the  Doppler  shift  of  the  atom,  then  after  each 
ott  tield  period  the  sign  ot  the  phase  error  with  respect  to 
the  e;  axis  will  be  reversed,  and  the  subsequent  *  pulse  will 
bring  the  phase  error  back  to  zero  The  phase  error  correc 
tion  cannot  be  made  exact  for  all  atoms,  of  course  because  of 
the  spread  in  the  precession  rate  about  e  caused  bv  the 
finite  velocity  spread  The  optimum  zero  field  period  tor 
minimizing  the  average  phase  error  for  that  half  of  the  distri 
button  which  is  being  cooled  is  roughly  r  i  A-./,  J1  where 
Awfi>  is  the  Doppler  width  of  the  transition 

This  phenomenon  is  similar  to  the  phase  spread  reduction 
that  occurs  in  photon  echoes  It  is  clear  that  for  a  given 
number  of  r  pulses,  the  net  accumulated  phase  error  and  the 
difference  in  the  accumulated  phase  error  lietween  one  side 
ot  the  distribution  and  the  other  will  be  respectively  less 
than  and  greater  than  their  values  for  the  case  of  no  zero 
fieid  period  Thus  the  use  of  the  zero  'ield  period  between  » 
pulses  will  permit  cooling  over  a  broader  velocity  distribu 
tion 

We  have  carried  out  a  theoretical  analysis  of  the  »  pulse 
cooling  technique  discussed  above  by  utibzing  the  quantum 
transport  equations  This  analysis  is  described  in  the  tol 
lowing  two  sections 

3.  THEORETICAL  FRAMEWORK:  THE 
QUANTUM  TRANSPORT  EQUATIONS 

For  the  purpose  of  completeness  we  dense  in  this  section 
the  tundamental  equations  describing  the  quantum  evolu 
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non  of  an  ensemble  of  two-level  atoms  interacting  w ith  resu 
nant  radiation  fields  The  equations  were  first  derived  hv 
Kolchenko  et  at  and  have  been  rediscovered  bv  several 
groups The  equations  include  the  effects  of  Doppler 
»hift.  grating  washout,  and  atomic  recoil  If  R  denotes  the 
center -of- mass  displacement  vector  and  r  the  internal  coor¬ 
dinates.  then  the  space-time  evolution  ot  an  atom  is  got 
erned  bx  the  Schrodinger  equation 
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where  the  tirst  term  on  the  right  hand  side  accounts  tor  the 
translational  motion  ot  the  atom  H  is  the  electronic  Hamu 
toman  and  V  i  R  r  1  is  the  interaction  potential  We  shall 
i  onsider  an  interaction  potential  tor  the  case  •  [  'he  induced 
elect r u  dipole  moment  of  an  atom  interacting  with  the  radi 
it  ion  tield 

W  e  begin  tn  expanding  R  r  •  ' as  a  linear  superposition 
•t  eigenstates  ♦...it  H 

f"R  r  ■  *  \  i  R  ■  ■  fexp-  -  -  . 

* here  the  prodabilitv  amplitude  i  .  R  jc  c  •unlstor’ne 
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Where  We  deptitv  R  »  R  ♦  R  and  R  »  R  -  R  . 
■Xppiv.r.c  'he  camr  procedure  le  using  expression  in 
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wnere  the  convective  derivative  p  T  V/  accounts  lor  the 
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need  not  be  swept,  and  the  axiling  becomes  exponential  with 
time  with  the  il  e>  cooling  time  given  by  Eq  <141.  The 
temperature  limit  most  commonly  assumed  for  the  steady - 
slate  laser  cooling  is  the  so-called  stochastic  limit  ikn T  » 
v  \  2* 

The  daahed  curves  in  Fig  .3  show  the  estimated  exponen¬ 
tial  cooling  tor  t  pulses  using  the  characteristic  cooling  time 
given  b\  expression  *  13i  The  net  cooling  time  is  seen  to  be 
dramaticalh  shorter  for  the  larger  initial  temperatures  but 
at  the  expense  of  considerable  higher  average  power  density 
In  this  connection  one  must  note  that  the  Lyman  . i  photons 
are  capable  ot  photoionmng  the  upper  state  of  the  transi¬ 
tion  The  rate  ol  energc  released  to  electrons  due  to  ioniza¬ 
tion  is  small  compared  w  ith  the  cooling  rates  tor  the  steady- 
state  cooling  power  levels  but  is  comparable  with  the  cooling 
rates  tor  the  r  pulse  average  power  levels  and  must  be 
watched  closelv  in  further  analysis  ot  the  practical  applica¬ 
tion  ot  -  pulse  cooling  Fhotoionization  of  the  upper  level 
also  gives  rise  to  phase  shifts  not  accounted  tor  in  the  two- 
level  treatment  ot  r  pulse  cooling  given  above,  hut  these 
-hitts  .an  be  expected  to  re  small  ow mg  to  the  much  smaller 
dipoie  matrix  element  between  the  upper  state  and  the  con¬ 
tinuum  .ompared  with  that  .«!  the  L.vman  o  transition 
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ABSTRACT 


A  t ime -dependent  theory  of  2-photon  induced  ioniza- 
n  is  presented.  Exact  analytical  expressions  are  obtained 
♦.he  total  population  of  hydrogenic-1  ike  systems  and  the 
ulation  of  the  intermediate  resonant  state  as  a  function 
•he  input  intensity  and  detuning  of  the  radiation  field, 
•eaporal  behavior  of  the  ionization  rate  is  calculated. 
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I.  INTRODUCTION 


Ifultiphoton  ionization  of  hydrogenic  systems  is 
concerned  with  the  dynamics  involved  in  the  coupling  of  external 
radiation  fields  with  bound  quantum  systems,  leading  to  the 
production  of  free  electrons  and  ions.  As  such,  the  subject  of 
multiphoton  ionization  has  been  the  subject  of  intense 
investigation  during  the  past  decade1.  One  of  the  simplest  but 
most  fundamental  problem  is  the  study  of  the  ionization  process 
when  the  input  radiation  field  induces  a  resonant  two-photon 
transition  into  the  continuum  from  the  ground  state  of  a 
hydrogenic-1 ike  atom.  In  this  case,  there  are  two  distinct 
regimes.  First,  if  the  input  radiation  is  not  in  resonance  with 
the  bound-bound  transition;  the  ionization  rate  can  be  computed 
by  using  perturbation  theory2 .  The  results  show  that  the  rate  for 
the  production  of  electron  is  proportional  to  the  square  of  the 
input  radiation  intensity.  However  in  the  case  in  which  the  input 
radiation  is  in  resonance  with  the  bound-bound  transition,  a 
simple  intensity  dependent  ionization  rate  can  not,  in  general, 
be  derived3 .  The  effect  of  saturation  of  the  bound  state 
population  difference  in  the  presence  of  the  resonance  radiation 
prevents  a  simple  interpretation  of  the  ionization  process.  The 
objective  of  this  work  is  to  present  a  self-consistent 
calculation  of  the  effects  of  resonance  radiation  on  the  two- 
photon  ionization  process.  The  approach  we  shall  take  is  to  use 


of  the  density  matrix  equations  in  the  computation  of  the 
temporal  evolution  of  the  atomic  density  in  the  presence  of 
radiation  fields.  Section  II  presents  a  detailed  discussion  and 
derivation  of  the  density  matrix  approach  which  takes  into 
account  the  ensemble  of  atoms  passing  through  radiation  field.  We 
show  in  Section  III  that  in  the  adiabatic  regime  ,  exact  time 
dependent  analytical  solutions  can  be  found  for  the  total 
population  of  the  atoms  as  well  as  for  the  ionization  rate.  The 
long  time  behavior  of  the  ionization  rate  shows  a  rather  simple, 
but  not  the  square  dependence,  of  the  input  intensity.  As  an 
example,  we  shall  apply  our  results  to  the  case  of  ionization  of 
a  H  beam  using  Lya  radiation.  And  finally  in  Section  IV  we 
summarize  the  main  results  of  our  work. 


II.  FORMULATION  OF  THE  PROBLEM 


The  starting  point  in  our  analysis  is  the  formulation  of 
a  set  of  density  matrix  equations  which  describes  the  model  shown 
in  Fig.  1.  The  bound  states  are  labelled  by  1 1  >  and  12  >,  while 
the  continuum  is  labelled  by  the  set  {  Ik  >  }.  The  resonance 
radiation  oscillates  at  frequency  0.  The  spontaneous  decay  rate 
from  the  state  1 2  >  to  the  ground  state  I  1  >  is  given  by  7.  In 
order  to  derive  a  set  of  equations  describing  the  evolution  of 
each  level  via  the  density  matrix  formalism,  it  is  important  to 
understand  the  behavior  of  the  coherent  coupling  between  state 
12  >  and  {  Ik  >  }.  This  is  accomplished,  in  a  rather 
straighf orward  manner,  by  studying  the  temporal  evolution  of  the 
wavefunction  in  state  12  >.  Consider  the  Schroedinger  equation 

d  l¥> 

ih  -  =  (  H0  -  n  E  )  l¥>  (1) 

d  t 

where  the  state  vector  !♦>  is  given  by 

!♦>  =  Cx  (t)  ll>  expC-iu/it)  +  C2  (t)  1 2>  exp(-iw2t) 

+  /  dk  C(k,t)  lk>  exp(-iwkt)  (2) 

such  that  1 1 > ,  1 2>  and  I k>  are  the  eigenstates  of  the  unperturbed 

Hamiltonian  He .  The  interaction  between  the  atom  and  the 


radiation  field  E  is  ~(i  E.  The  near  resonant  coupling  between 
state  I 2>  and  the  continuum  can  be  obtained  by  projecting  the 
state  vector  !♦>  into  the  respective  space.  They  are 

d  C2(t) 

ih  -  =  -  /i2i  B  exp(iw21t) 

d  t 

-  /  dk  /*2k  E  C(k,t)  exp(-iwk2t)  (3a) 


d  C(k,t) 

ih  -  =  -  fik2  E  C2(t)  exp(iwk2t)  (3b) 

d  t 

where  jiap  and  uap  are  the  electric  dipole  matrix  element  and 
frequency  separation  between  states  I a>  and  l/?>;  respectively.  In 
the  derivation  of  Eqs . (3) ,  we  assume  that  the  coupling  among  the 
continuum  states  is  negligible  and  the  off -resonant  two-photon 
transition  between  the  ground  state  and  the  continuum  has  a  small 
contribution  in  comparison  with  the  resonant  transition. 

Equations  (3)  can  be  reduced  further  by  formally 
integrating  Eq. (3b)  and  use  the  result  in  Eq.(3a).  This  procedure 
leads  to  the  following  expression  for  the  evolution  of  C2 (t) 


d  C2(t) 


ih -  =  -  /*2lE  Cx(t)  exp(iw21t) 

d  t 

+  (ih)-1  ;  dk  /  dt*  M2 i,E(t)/ih2E(t’) 

x  exp(-iwh2 C2 (t ’ )  (4) 

Now  consider  the  approximations  involved  in  the  further 
simplification  of  Eq. (4) .  First,  we  shall  assume  that  the 
amplitude  C2(t*)  is  slowly  varying  on  the  time  scale  of  the 
radiation  intensity  and  1  /a/  k  2  •  Then  we  can  extract  C2(t’)  out  of 
the  integral  and  replace  it  by  C2 (t) .  Hence  the  integral  involves 
only  the  intensity  of  the  radiation  field  and  the  medium 
parameter.  This  integral  is  complex  and  the  expression  can  be 
written  as  the  sum  of  a  real  and  an  imaginary  part;i.e. 

(ih)'1  /  dk  /dt’  pak*  E(t)Mk2«E(t')  exp  (-iwk  2 (t-t ’ ) ) =-ih  (I*-iA)  . 
Hence  the  net  effect  of  the  coupling  betwen  state  I  2>  and  the 
continuum,  in  the  adiabatic  limit,  is  the  appearance  of  an 
irreversible  transfer  of  energy  (  intensity  dependent  )  from  the 
bound  state  to  the  continuum.  Furthermore,  there  exists  an 
intensity  dependent  level  shift  of  I2>.  We  shall  write  the 
intensity  dependence  of  T  and  A  as  (7 1  and  61;  respectively.  This 
result  provides  the  starting  point  in  the  derivation  of  the 
density  matrix  equations.  However  if  the  assumption  that  the 
amplitude  C2  is  not  slowly  varying,  then  additional  complications 
may  arise  which  may  lead  to  the  existence  of  two-photon  Rabi 
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flopping,  sven  in  the  presence  of  the  continuum.  This  fact  was 
pointed  out  by  Kumenov  and  Perel4  in  their  discussion  of  single¬ 
photon  induced  ionization.  We  have  carried  out  similar 
calculations  for  the  case  of  two-photon  resonant  ionization  and 
we  found  that  the  ionization  process  can  be  eliminated  for 
certain  range  of  the  radiation  intensity. 

Let  p n ,  p 22  and  P\  2  be  the  population  of  state  ll>, 
population  of  state  I 2>  and  the  optical  coherence  between  states 
I 1>  and  I 2> ;  respectively.  The  equations  of  motion  is  given  by 
d  P 11 

ih  =  7P22  ~  Mis  E  psi  +  p12  P21  E  (5a) 

d  t 

d  P 22 

ih  _  ~^p22  -  a  I  p  22  -  Mu  ®  Pis  +  Psi  Mis  ®  (5b) 

d  t 

d  Pis  7 

ih  =  -hw2j  P\ 2  _ih  P12  -  Mis  E  (  pz 2  ~  P11  )  (5c) 

d  t  2 

where  w21  =  u9  +  61.  The  set  of  equations  (5)  describes  the 
evolution  of  the  bound  systems  in  the  presence  of  radiation 
fields.  It  is  interesting  to  note  that  Eqs . (5)  describe  the 
dynamics  of  the  model  shown  in  Fig. 2.  Hence  within  the 
approximations  stated  above,  the  dynamics  of  two-photon 


ionization  process  can  be  described  in  terms  of  the  dynamics  of 
am  open  two-level  atom,i.e.  the  atom  is  coupled  with  the 
continuum  which  is  treated  as  a  reservoir.  In  the  next  section  we 
shall  solve  these  equations  and  obtain  the  probability  that  the 
atoms  are  ionized  by  means  of  the  two-photon  excitation. 


III.  SOLUTIONS 


The  set  of  density  matrix  equations  have  exact 
analytical  solutions  provided  that  the  following  assumptions  are 
made  concerning  the  temporal  evolution  of  the  radiation  field, 
the  population  difference,  p2 2  -  p11(  and  the  optical  coherence. 
First,  the  radiation  field  is  assumed  to  have  a  pulse  duration 
that  is  longer  than  any  response  time  of  the  medium.  For  sake  of 
understanding  of  the  process  of  resonant  two-photon  ionization, 
the  radiation  field  is  assumed  to  be  cw.  Second,  the  population 
difference  is  assumed  to  be  slowly  varying  in  a  time  compared  to 
the  smaller  of  the  natural  linewidth  or  laser  detuning  from 
intermediate  resonance.  This  is  just  the  Rate  Equation 
Approximation.  And  last,  the  laser  detuning  from  resonance  is 
assumed  to  be  small  enough  so  that  one  can  neglect  the  Bloch- 
Siegert  shifts.  This  last  assumption  entails  making  the  Rotating 
Wave  Approximation. 

Taking  these  assumptions  into  account,  let  us  write  the 
electric  field  as  E  =  E0  cos  At  and  the  on-resonance  Rabi 
flopping  frequency  as  R0  =  p12  Ee/2h.  In  the  rotating  wave 
approximation,  the  optical  coherence  p12  is  given  by  p 12=  p12 
exp(ifit) .  Then  the  density  matrix  equations  become 


<*  p  11 

=  7  P22  +  i  R*  P 21  “  i  Re*  P 12  (6a) 

d  t 

d  p22 

=  -7  p22  -  a  I  p22  +  i  R0*  p12  -  i  R,  p21  (6b) 

d  t 

d  p12 

-  +  [iA  +  7/2]  p12  =  i  Re  (  P22  -  Pxi  )  (6c) 

d  t 

Integration  of  Eq. (6c)  and  using  the  Rate  Equation  and 
Rotating  Wave  Approximations  yield  the  solution  for  p12.  The 
expression  for  p21  can  be  obtained  by  taking  the  complex 
conjugate  of  px 2.  They  are  given  by 

P12  =  i  Ra  L12  (  p22  -  Pi  1  )  (7a) 

P 21  =  “i  Re*  kl2*  (  p22  ~  P11  )  (7b) 
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where  L12  =  {  iA  +  7/2  } * 1  is  the  complex  Lorenztian  lineshape 
f  actor . 


Using  Eqs .  (7a)  and  (7b)  in  Eqs .  (6a)  and  (6b)  yields  a 
close  system  of  equations  for  plx  and  p2 2. 

d  Pi  x 

=  7  P 22  +|R*I2  (P 22  ~  P11)  (  ^*12  )  (8a) 

d  t 


d  p22 

=  -(7  +  a  I  )p22  -  IRg  I2  (p  2  2  -  P11)  (L12)  (8b) 

d  t 


where  we  shall  assume  that  the  initial  conditions  are  such  that 
Pil  (  )  =  Nfl  and  p22  (  t=-«  )  =  0.  N0  is  the  initial  density 

of  the  atoms  in  the  ground  state.  Im  (LI2)  standa  for  the 
imaginary  part  of  L12. 

The  reduced  set  of  equations  (8)  has  the  following  exact 


analytical  solutions 


N#  K 

Pjj  =  -  exp>  5(7  *  °  1  '  K'tj  »mh(  X,t  )  (©a) 

X  j 

Pii  =  -  P77  O  I  /*  ^t  P?7  (^  )  (8b) 

where  K  =|R.|2  Im  (L12)  and  Xj={(7  ♦  oI+  K) 2  -  4cIK}1/2 

Eqs . (9)  constitute  one  of  the  two  main  results  of  this 
work  and  they  are  valid  in  the  regime  where  the  three  fundamental 
assumptions  are  satisfied.  Figure  3  shows  a  plot  of  the  temporal 
evolution  of  the  excited  state  population  density.  It  shows  that 
the  population  of  the  excited  state  achieves  a  large  value  within 
a  short  period  of  time  due  to  the  Rabi  flopping  process.  After 
that. the  single  photon  ionization  process  takes  place  to  deplete 
the  excited  state  population.  A  nearly  complete  depletion  occurs 
within  few  hundreds  of  the  spontaneous  emission  time.  Figure  4 
shows  the  temporal  behavior  of  the  total  population  density  and 
exhibits  the  predicted  depletion  due  to  the  ionization  process  in 
the  high  intensity  regime;  i.e.  500  Kw/cm2.  However,  if  the 
intensity  is  restricted  to  1  Kw/cm2,  negligible  ionization  takes 
place.  Hence  the  integrity  of  the  atomic  beam  is  preserved. 

Figure  5  shows  the  behavior  of  the  total  population  density  for 
fixed  intensity  (  500  Kw/cm2  )  as  a  function  of  detuning  from 
line  center.  For  large  detuning  the  ionization  process  becomes 
small.  This  behavior  is  understood  in  terms  of  the  nonresonant 


nature  of  the  detuned  laser  action.  These  plots  consider  the  case 
of  hydrogen  atoms  that  are  being  ionized  by  the  action  of  a 
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and  2p  states  of  hydrogen  The  cross  »e  • : or  for  ionization  from 

2p  to  the  continuum  was  assumed  to  he  of  *he  c'rder  of  10  ‘f  cm* 
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A  quantity  of  interest  is  the  resonant  two  photon 
ionization  rate  It  was  stressed  that  this  quantity  has  an 
intuitive  physical  meaning  provided  that  perturbation  theory 
holds,  1  e  the  Fermi  s  golden  rule  is  used  in  the  discussion  If 
we  were  to  approach  such  a  quantity  in  the  regime  of  saturation 
of  the  bound  transition,  then  the  ionization  rate  is  defined  as 


d  (  Pi  1  4  P?  ?  ) 


P 11  *  P 77  dt 

Using  the  exact  analytical  expressions  (9) , 
photon  ionization  rate  is  given  by 


(10) 


the  resonant  two- 


a  I  K  1 

R.on  = -  -  (ID 

X  2  coth(\2t)  +  (7  +  (7 1  +  K)/2\2 

It  is  interesting  to  note  that  the  rate  is  time  dependent  due  to 
the  fact  that  the  exact  solutions  describe  the  temporal  behavior 
of  the  population  density  as  a  function  of  the  laser  intensity. 
The  long  time  behavior,  i.e.  as  t  tends  to  infinity,  is  given  by 


2  al  K 


(12) 


Rion  ~ 

2X2  +  7  +  al  +  K 

which  is  a  nonlinear  function  of  the  laser  intensity,  contrary  to 
the  results  obtained  via  the  perturbation  regime.  Eq. (12)  is  the 
second  of  the  two  main  results  of  this  work.  It  shows  that  the 
two-photon  ionization  rate  reduces  to  the  perturbation  results 
in  the  limit  of  low  intensity  and  zero  detuning  from  the 
intermediate  state  resonance , i . e .  the  ionization  rate  is 
proportional  to  I2.  While  in  the  high  intensity  limit  the  rate 
becomes  linear  in  intensity. 


.J  -  4  , 


t j  u  i.i  *-4  »_>  t.t  »-i  «a  >j  i.<  fc-A  >i  ».t 


IV.  SUMMARY 


We  presented  a  theory  of  resonant  two-photon  ionization 
of  hydrogenic-like  systems.  The  analysis  shows  that  exact 
analytical  results  are  obtained  for  the  total  population  of  the 
atoms  and  the  ionization  rate,  provided  that  three  fundamental 
approximations  are  satisfied.  They  are  the  slowly  varying 
envelope,  rate  equation  and  rotating  wave  approximations.  These 
approximations  are  valid  for  situations  where  the  laser  pulse 
duration  is  long  compared  to  the  characteristic  times  of  the  atom 
and  that  the  laser  is  tuned  near  resonance  to  the  bound-bound 
transition . 

The  ionization  rate  was  found  to  be  a  nonlinear  function 
of  the  laser  intensity  due  to  the  saturation  process  of  the  bound 
state  transition.  It  reduces  to  the  result  obtained  from 
perturbation  theory  in  the  limit  of  low  intensity.  In  the  high 
intensity  limit  the  rate  becomes  linear  in  the  intensity,  which 
reflects  the  fact  that  the  saturation  process  plays  an  important 
role  in  determining  the  dynamics  of  the  two-photon  ionization 


process 


In  the  course  of  this  work,  we  found  that  under  certain 


specific  conditions  it  is  possible  to  change  the  continuum  in 
such  a  manner  that  ionization  does  not  take  place.  Preliminary 
work  points  to  the  criterion  that  the  laser  pulse  must  be  short 
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enough  that  the  adiabatic  condition  is  no  longer  valid.  In  this 
case  it  is  not  possible  to  extract  C2  from  the  integral  in  Eq. 

(4) .  It  can  be  shown  that  there  exists  a  range  of  parameters  ( 
laser  pulse  duration  and  intensity  )  such  that  state  |2>  does  not 
experience  an  irreversible  energy  transfer.  This  results  implies 
that  the  radiation  field  has  modified  the  continuum  such  that  the 
ionization  process  does  not  take  place.  We  shall  expose  this 
subject  matter  in  full  detail  in  another  publication. 

And  last,  a  brief  comment  is  in  order  concerning  the 
role  of  ionization  in  the  proposed  laser  cooling  of  a  H  beam. 
Assuming  that  it  takes  of  the  order  of  microsecond  to  achieve  a 
divergence  of  the  order  of  /irad  and  the  Lya  laser  is  saturating 
the  Is  to  2p  transition  (  for  optimum  radiation  cooling  ) ,  then 
Fig.  2  shows  that  a  significant  fraction  of  atoms  are  ionized 
during  the  cooling  process.  This  calculation  uses  a  single-photon 
cross  section  of  10" 18  cm2  for  the  transition  from  the  2p  to  the 
continuum.  A  more  careful  analysis  of  the  cross  section  must  be 
done  in  order  to  assess  correctly  the  role  of  two-photon  induced 
ionization  on  the  laser  cooling  process. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


.  The  model  used  to  discuss  the  dynamics  of  two-photon 
resonant  ionization.  0  is  the  frequency  of  the  input 
radiation  field  and  7  is  the  spontaneous  emission  decay 
rate . 

.  The  reduced  model  showing  that  the  net  effect  of  the 
interaction  between  state  |2>  and  the  continuous  is  to 
introduce  an  incoherent  decay  rate  a  I  and  a  shift  in 
the  energy  level  of  the  intermediate  state. 

.  Temporal  evolution  of  the  excited  state  population 
in  the  presence  of  two-photon  induced  ionization. 

.  Temporal  evolution  of  the  total  density  of  atoms  in 
the  presence  of  two-photon  induced  ionization 
The  Lya  laser  is  tuned  on-resonance. 

.  Temporal  evolution  of  the  total  density  of  atoms  in 
the  presence  of  two-photon  ionization,  for  fixed 
laser  intensity.  The  Lya  radiation  intensity  is  set 
at  500  kW/cm2 . 
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APPENDIX  C 

Velocity-specific  atomic-state  selection  in  an  atomic  beam  by 
continuous-wave  optical  pumping 
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We  preaenl  experimental  reeulta  demonstrating  strong  velocity-distribution  modification  in  e  specific  atomic  state 
using  optical  pumping  between  hyperfine  split  levels.  TTiess  results  show  that  the  velocity  of  a  specific  state  can 
be  defined  within  6%.  corresponding  to  a  local  temperature  of  1  K.  The  velocity-epeeifie-etate  selection  is  accom¬ 
plished  by  using  a  nearly  effusive  atomic  beam  of  sodium  and  a  single  frequency-stabilised  dye  laser. 


In  this  Letter  we  describe  experimental  measure¬ 
ments  demonstrating  the  use  of  optical  pumping  for 
-generating  an  atomic  beam  with  one  of  the  ground  states 
of  the  atomic  species  characterized  by  a  nearly  mono¬ 
chromatic  longitudinal  velocity  distribution.  Such  a 
technique  eliminates  the  complications  associated  with 
using  a  rotating -drum  velocity  selector.1  The  results 
show  the  significance  of  velocity -specific  optical 
pumping  and  velocity-distribution  modification,  which 
characterize  many  nonlinear  spectroscopy  experi¬ 
ments.  ^ 

The  proposed  basic  approach  was  discussed  earlier.4 
We  have  implemented  a  slight  modification  to  improve 
the  original  proposal,  and,  for  the  sake  of  completeness, 
we  review  the  basic  principles  here.  The  technique  is 
based  on  the  optical-pumping  properties  of  atomic 
states  split  by  hyperfine  interaction.  In  these  experi¬ 
ments,  we  used  the  D2  line  of  atomic  sodium  at  589  nm. 
An  energy-level  diagram  of  sodium  is  shown  in  Fig.  1, 
in  which  we  have  included  only  the  3 s  2Sm~3p  2Py 2 
transition.  The  total  angular-momentum  quantum 
number  F  is  the  sum  of  the  nuclear  spin  (/  -  3/2  for 
*'Na)  and  the  particular  J  value  for  the  given  level. 
The  various  hyperfine  splittings  are  shown  in  Fig.  1. 
For  these  transitions,  the  relevant  dipole-selection  rules 
are  AF  «  0,  ±1.  Hence  we  see  that  an  atom  in  the  3p 
2Fn/2<F  ”  3)  state  can  decay  by  spontaneous  emission 
only  to  the  3«  2S m (F  ”  2)  state.  Similarly,  an  atom  in 
the  3p  2Pm(F  m  0)  state  can  decay  by  spontaneous 
emission  only  to  the  3s  2S  m(F  m  1)  state.  However,  the 
3p  2Pyi(F  •  I)  state  and  the  3 p  2Pyi(F  *  2)  state  can 
decay  spontaneously  to  either  of  the  3s  2S\n  ground 
states.  Therefore  an  atom  initially  in  the  F  -  2  ground 
state  that  is  excited  to  the  F  m  2  or  F  ■  1  upper  state  can 
end  up  in  the  F  -  1  ground  state.  Hence,  by  optical 
pumping,  the  entire  atomic  population  can  be  trans¬ 
ferred  from  the  F  *  2  ground  state  to  the  F  -  1  ground 
state.  In  Fig.  1  the  transitions  that  are  not  amenable 
to  optical  pumping  are  shown  with  solid  lines,  and  the 
dotted  lines  show  transitions  that  can  be  optically 
pumped.  By  an  appropriate  arrangement  of  optical 
frequencies  and  geometry,  it  is  possible  to  convert  the 
original  broad  velocity  distribution  of  one  of  the  ground 
states  to  a  narrow  velocity  distribution.  To  see  how  this 


is  accomplished,  we  consider  the  optical-pumping 
scheme  shown  in  Fig.  2.  The  energy-level  spacing  in  the 
ground  state  is  sufficiently  small  that  one  can  assume 
that  both  ground  states  are  thermally  populated.  We 
consider  an  optical-beam  arrangement  in  which  the 
atomic  beam  is  sufficiently  collimated  that  we  need  not 
consider  any  transverse  velocity  effects.  We  further 
assume  that  the  initial  longitudinal  velocity  distribution 
is  thermal.  The  optical-pumping-induced  velocity- 
distribution  modification  is  a  two-step  process.  The 
first  step  requires  the  frequency  of  a  laser  beam  prop¬ 
agating  perpendicular  to  the  atomic  beam  (to  avoid 
Doppler  effects)  to  be  resonant  with  the  3s  2S\n(F  * 
2)-3p  2Pin(F  *  2)  [or  3p  2Pm(F  •  I)]  transition  (e.g., 
transition  A  of  Fig.  2).  Assuming  that  the  laser  inten¬ 
sity  is  comparable  with  the  transition-saturation  in¬ 
tensity  (given  by  hu/oT  1,  where  a  is  the  absorption 
cross  section  and  T ,  is  the  upper-state  lifetime)  and  that 
the  transit  time  is  long  compared  with  T 1,  the  atoms  in 
the  F  =  2  ground  state  will  decay  eventually  to  the  F  * 

1  ground  state,  where  they  will  no  longer  be  resonant 
with  the  laser  light  (being  out  of  resonance  by  some  1700 
MHz).  The  second  step  in  the  velocity-distribution 
modification  process  is  to  pick  a  velocity  class  (i.e.,  a 
group  of  atoms  moving  within  some  velocity  Ac  about 
u )  in  the  F  ■  1  ground  state  and  transfer  it  to  the  F  * 

2  ground  state.  This  is  accomplished  by  taking  ad- 
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Fig.  1.  Energy-level  diagram  for  the  0\  line  of  atomic  sodium 
at  589  nm  showing  the  hyperfine  splitting. 
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Fig.  2.  Left,  atomic-energy- level  scheme  showing  the  opti¬ 
cal-pumping  sequence  for  velocity  selection.  Right,  geometry 
for  the  optical-pumping  laser  beam  and  the  velocity-selection 
laser  beam. 


vantage  of  the  longitudinal  Doppler  shift  A  beam  of 
light  at  frequency  u>  is  directed  at  some  angle  4>  with 
respect  to  the  atomic  beam.  In  this  case,  we  want  to 
choose  a  specific  velocity  group  such  that  it  is  in  reso¬ 
nance  with  the  laser  on  the  3s  2Sui(F  m  1)— 3p  2P$n(F 
=  1)  (or  3 p  2P^(F  *  2)]  transition  at  wo'  (or  wo')-  The 
resonance  condition  is  given  by  w  -  wo  -  K  •  v  *  0, 
where  K  is  the  K  vector  of  the  laser  beam,  wo  *  w0'  or 
wo",  and  u  is  the  atomic  velocity  of  the  velocity  group  we 
wish  to  excite.  When  w  satisfies  this  resonance  condi¬ 
tion,  all  atoms  moving  in  the  velocity  group  u  will  be 
excited.  Transition  B  of  Fig.  2  shows  an  example  when 
wo  *  wo'.  Assuming  again  that  the  laser  intensity  is 
comparable  with  the  transition-saturation  intensity,  we 
find  that  those  atoms  excited  will  eventually  decay  back 
to  the  F  ■  2  ground  state,  where  they  will  be  out  of 
resonance  again  by  some  1700  MHz.  At  this  point,  the 
velocity  distribution  of  the  F  *  2  ground  state  has  been 
modified  from  a  broad  thermal  distribution  to  a  narrow 
distribution  moving  at  velocity  v.  An  example  of  the 
geometry  is  shown  in  Fig.  2. 

The  width  of  the  velocity  distribution  is  limited  in 
principle  by  the  natural  linewidth  of  the  transition,  i.e., 
the  beam  will  pump  all  atoms  whose  velocity  for  a  given 
w  and  wo  satisfies  the  inequality  jw  —  wo  -  K  ■  t/|  <  y, 
where  2y  is  the  transition  linewidth  (FWHM  in  ra¬ 
dians).  In  this  particular  example,  this  limit  cannot  be 
realized  because  there  are  two  optical-pumping  tran¬ 
sitions  involving  both  the  F  *  1  and  F  m  2  levels  in  the 
upper  state,  yielding,  in  principle,  a  velocity  doublet  If 
their  frequency  separation  is  given  by  <5,  then  the  best 
longitudinal  monochromaticity  that  can  be  realized  is 
estimated  by  replacing  y  with  i  in  the  above  inequality, 
showing  that  the  beam  monochromaticity  should  be  of 
the  order  Ao/o  ~  i/(w  -  wo),  where  6  is  in  radians.  For 
these  experiments  in  sodium  we  estimate  that  the 
minimum  £w/v  is  of  the  order  of  3%  (for  w  tuned  to  the 
3s  2Si/a(F  ■  2)-3p  2Pvs(F  ■  3)  transition]. 

The  experimental  configuration  to  demonstrate  this 
behavior  is  shown  in  Fig.  3.  The  two-step  process  is 
accomplished  with  a  single  stabilized  tunable  dye  laser 
tuned  to  539  nm.  This  technique  was  developed  not 
only  to  enable  a  single  dye  laser  to  accomplish  the  ve¬ 
locity-distribution  modification  but  also  to  allow  that 


dye  laser  to  be  used  as  the  radiation  source  for  future 
experiments  on  the  prepared  monochromatic  atomic 
beam.  For  those  experiments,  the  dye  laser  is  operated 
in  resonance  with  the  3a  2Sm(F  m  2)-3p  2P$n(F  ■  3) 
transition.  Hence  these  experiments  were  performed 
with  the  laser  tuned  to  that  resonance.  The  frequency 
offset  for  adjusting  the  laser  frequency  to  be  in  reso¬ 
nance  with  the  3a  2S\n(F  ■  2)-3p  2Py2(F  ■  2)  transi¬ 
tion  was  provided  by  an  acousto-optic  modulator  tuned 
to  59.6  MHz.  The  downshifted  optical  beam  was  di¬ 
rected  into  the  atomic  beam  at  90s  to  the  direction  of 
atomic-beam  propagation.  This  beam  transferred 
nearly  all  the  population  of  the  F  *  2  ground  state  into 
the  F  ■  1  ground  state.  The  second  step  was  accom¬ 
plished  by  taking  the  original  laser  beam  [tuned  to  the 
3s  2Sin(F  ■  2)-3p  2Pyi(F  »  3)  resonance]  and  directing 
it  to  intersect  the  atomic  beam  at  45°.  Atoms  moving 
at  roughly  1.4  X  10s  cm/sec  will  see  this  beam  as  being 
on  resonance  with  the  intermediate  transition  given  by 
3s  2Sin(F  «  1) — 3s  2Sin(F  «  1)  [or  3p  2P3/2(F  -  2)]. 
These  atoms  will  be  pumped  into  the  F  ■  2  ground 
state.  In  this  manner  the  atoms  in  the  F  *  2  ground 
state  are  all  moving  at  nearly  the  same  longitudinal 
speed. 

To  verify  the  resultant  velocity  distribution  experi¬ 
mentally,  a  second  frequency-stabilized  tunable  dye 
laser  was  directed  along  the  axis  counterpropagating  to 
the  atomic  beam,  as  shown  in  Fig.  3.  A  photomultiplier 
detected  the  scattered  fluorescence  from  the  second 
laser  at  a  position  downstream  from  the  optical- 
pumping  and  velocity-selection  region.  In  the  absence 
of  any  optical  pumping  or  velocity  selection,  we  expect 
a  double-humped  fluorescence  spectrum  as  the  laser  is 
tuned  through  the  resonances.  The  data  are  shown  in 
Fig.  4(a).  The  first  peak  is  the  Doppler- broadened 
spectrum  corresponding  to  all  transitions  originating 
out  of  the  F  *  2  ground  state,  and  the  second  peak 
corresponds  to  the  F  *  1  ground  state.  These  peaks  are 
separated  by  approximately  1700  MHz.  Figure  4(b) 
shows  the  effect  on  the  fluorescence  spectrum  as  the 
optical-pumping  beam  is  turned  on.  The  peak  corre¬ 
sponding  to  the  F  m  2  transitions  is  reduced  while  the 
peak  corresponding  to  the  F  *  1  transitions  is  increased. 
Because  of  the  geometry  of  the  probing  beam  and  the 
fact  that  the  probe  beam  was  larger  than  the  optical- 
pumping  beam,  not  all  atoms  that  are  probed  have  been 
pumped.  Hence  the  F  •  2  level  does  not  appear  to  be 
completely  depleted  in  these  data.  However,  a 
transverse  probe  beam  (not  shown  in  Fig.  4)  is  capable 
of  measuring  the  remaining  population  (with  the  loss 
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Fig.  3.  Experimental  configuration  for  optical  pumping  and 
velocity  selection. 
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Fig.  4.  Probe  scan  us'ng  the  second  dye  laser  to  monitor  the 
Doppler-broadened  fluorescence  spectrum,  (a)  The  spectrum 
obtained  with  no  optical-pumping  beam  or  velocity-selection 
beam  present.  The  double  hump,  separated  by  roughly  1700 
MHz.  is  due  to  the  hyperfine  splitting  in  the  ground  state,  (b) 
In  the  presence  of  optical-pumping  beam  A  (see  Fig.  2).  we 
observe  a  depletion  of  the  F  »  2  ground  state  and  an  en¬ 
hancement  of  the  F  *  l  state.  Ic)  Fluorescence  spectrum  in 
the  presence  of  both  the  optical -pumping  beam  and  the  ve¬ 
locity-selection  beam.  The  spike  on  the  low-frequency  side 
shows  the  clear  evidence  of  state-specific  velocity  selection. 
A  high-resolution  display  of  this  peak  shows  a  width  of  90 
MHz,  corresponding  to  a  velocity  homogeneity  of  5%. 


of  velocity  resolution)  and  shows  nearly  complete 
transfer  of  population  from  level  2  to  level  1.  Figure 
4(c)  shows  the  effect  on  the  velocity  distribution  when 
the  velocity-selection  beam  is  turned  on  at  $  *  45®  with 
respect  to  the  oncoming  atomic  beam.  A  strong  peak 
is  observed  corresponding  to  atoms  in  the  1.4  X  10s- 
cm/sec  velocity  group.  A  high-resolution  display  of  this 
peak  shows  a  width  of  90  MHz,  implying  a  velocity  ho¬ 
mogeneity  of  6%,  in  good  agreement  with  the  estimates 
given  above.  The  equivalent  local  temperature  is  of  the 
order  of  1  K.  The  velocity  homogeneity  was  limited  by 
power  broadening  and  diffractive  divergence  of  the 
velocity-selecting  optical  beam,  thus  eliminating  the 
possibility  of  observing  the  velocity  doublet. 

The  amplitude  of  the  velocity  peak  is  larger  than 
expected  based  on  the  data  of  Fig.  4(a)  because  of  op¬ 
tical  pumping  effects  caused  by  the  weak  probe  beam 
in  collinear  geometry.4  When  a  transverse  probe  is 
used,  the  increased  signal  amplitude  is  in  agreement 
with  that  expected  from  the  transferred  population  that 
is  due  to  the  velocity-selecting  beam. 


In  these  experiments  we  have  demonstrated  the  ef¬ 
fects  of  velocity-specific  atomic-state  selection  by  op¬ 
tical  pumping.  By  incorporating  an  acousto-optic 
modulator,  we  have  been  able  to  accomplish  the  velocity 
selection  while  the  laser  is  tuned  to  a  nonoptically 
pumped  transition.  Hence  additional  experiments  may 
be  performed  by  using  this  atomic  beam  and  the  same 
laser  with  the  laser  directed  toward  the  atomic  beam  at 
90*.  This  eliminates  the  problem  of  using  the  laser 
tuned  to  an  intermediate  transition  and  then  per¬ 
forming  the  experiments  on  the  pumped  transition  by 
adjusting  the  angle  of  the  beam  away  from  90* ,  as  sug¬ 
gested  by  Stroud.4  This  is  important  for  experiments 
in  which  standing-wave  effects  occur,  such  as  in  the 
optical  Stern-Gerlach  effect,8-7  which  we  are  now  ex¬ 
amining.  The  same  laser  can  now  be  used  to  replace  the 
classic  hot-wire  detector  used  to  measure  atomic-beam 
profiles.  Since  this  laser  is  tuned  to  the  3s  2S\n(F  * 
2>— 3 p  2Pm(F  *  3)  transition,  it  can  be  used  as  the 
atomic-beam  detector  downstream.  In  our  experi¬ 
ments.  the  laser  is  focused  to  less  than  30  nm  by  using 
a  cylindrical  lens  and  is  mechanically  scanned  to  provide 
the  transverse  profile  of  the  atomic  beam.  Because  of 
the  frequency  resolution,  it  samples  only  atoms  in  the 
F  *  2  state.  A  photomultiplier  is  then  used  to  detect 
the  side-scattered  fluorescence.  In  this  manner  high¬ 
speed  atomic- beam  detection  is  accomplished  with 
spatial  resolution  comparable  with  that  obtained  using 
a  hot-wire  detector  with  considerably  increased  band¬ 
width  and  none  of  the  difficult  technological  problems 
associated  with  hot  wires. 
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H~  Photoneutral ization 


Photoneutralization  presents  an  attractive  alternative  for 
neutralization  of  H~  beams  because  of  the  possibility  of 
drastically  decreasing  the  divergence  of  the  resultant  neutral 
beam.  Advances  in  laser  technology  make  photoneutralization 
increasingly  attractive  from  the  point  of  view  of  space  and  power 
requirements  as  well.  One  of  the  traditional  difficulties 
associated  with  weapons  use  of  photoneutralization  is  the 
process’s  low  cross  section.  This  presentation  presents  some 
ideas  for  overcoming  this  limitation.  Assessment  of  the 
technical  feasibility  of  these  techniques  require  experimental 
studies  which  should  be  initiated  shortly. 

The  first  possibility  I  will  address  for  enhancing 
photoneutralization  cross  sections  is  to  take  advantage  of  two 
resonant  structures  in  the  cross  section  near  the  n=2  threshold 
of  atomic  hydrogen,  a  reshbach  and  a  shape  resonance.  The 
Feshbach  resonance  is  located  just  below  the  n=2  threshold  and 
has  a  cross  section  which  is  predicted  to  be  about  300  times  the 
background  cross  section  at  that  energy.  While  it  has  been  seen 
experimentally,  previous  experiments  have  insufficient  resolution 
to  determine  details  about  its  height  and  width. 

The  second  resonance  of  interest  is  the  shape  resonance  which 
about  31  meV  above  the  n*2  threshold.  Although  its  cross 
section  is  predicted  to  be  below  that  of  the  Feshbach  resonance 
by  about  an  order  of  magnitude,  calculations  indicate  that  this 
reaction  strongly  favors  leaving  the  residual  hydrogen  atom  in 
the  n*2  state,  giving  the  ejected  electron  a  recoil  energy  of 
only  31  meV.  Beam  sensing  is  then  possible  at  visible 
wavelengths  via  the  2p  to  n*3  transitions.  Additionally, 
divergence  due  to  the  recoiling  electron  is  extremely  small. 
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A  third  technique  which  should  be  investigated  is  the 
possibility  of  electric  and  magnetic  field  induced  resonances. 
When  the  hydrogen  ions  are  placed  in  an  external  field, 
metastable  excited  states  above  the  ionization  threshold  become 
possible.  Since  the  transition  becomes  bound-bound  rather  than 
bound-free,  the  wave  function  overlap  (and  hence  the  cross 
section)  is  expected  to  increase  substantially.  Using  an 
external  magnetic  field  has  a  possible  drawback  of  the  fringing 
fields  introducing  a  twisting  motion  in  the  beam  which  may  add  to 
beam  divergence.  This  could  be  circumvented  by  using  an  electric 
field  instead.  Detailed  study  of  electric  field  induced 
resonances  needs  to  be  done  before  more  can  be  said  about  the 
details  of  this  enhancement. 

In  conclusion,  photoneutralization  of  H  beams  promises  to  be 
an  attractive  long-term  alternative  to  other  NPB  neutralization 
techniques.  This  process  presents  the  possibilities  of  lower 
contributions  to  beam  divergence.  The  disadvantage  of  low  cross 
sections  which  require  large  laser  powers  may  be  offset  by 
successful  utilization  of  resonant  schemes  which  can  enhance  the 
cross  section  substantially.  Detailed  investigations  of  these 
processes  should  begin  immediately  so  that  the  technology  can  be 
tested,  understood,  and  developed  if  feasible. 
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OUTLINE 


(1)  RL  Facilities 

— High  Power  Narrow  badn  tunable  VUV 
--Hydrogen  Beam  Facility 

(2)  Resonant  Photoneutralisation 
— Feshbach  Resonance 
--Shape  Resonance 

(3)  Field  Enhanced  Photoneutralization 


(4)  Summary 


TUNABLE  TRANSFORM  LIMITED 
BANDWIDTH  VUV  SOURCE 

\  =  1 216& 


FESHBACH  RESONANCE 


--Relatively  Large  Cros*  Section 

--Is  Final  State 

10  eV  e ~  recoil  produces  5/irad  divergence  in  200  MeV 
beam 

--Previous  Work 

electron  scattering 
LAMPF  experiments 

— Future  Work 

High  Resolution  studies 
Resolve  other  states 


D-8 


SHAPE  RESONANCE 


-smaller  cross  section 

-lower  divergence 

0-3  ad  with  a  200  MeV  H~ 

■previous  work 

Bryant,  et  al .  (LAMPF) 
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SUMMARY 


— Feshbach  Resonance 


High  cross  section  /  highest  divergence 


— Shape  Resonance 


Lower  cross  section  /  lower  divergence 


--Field  Enhanced  Photoneutralization 


NOTE:  studies  can  be  done  now  on  a  benchtop  with  existing 


laboratory  equipment! 


